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JUGFAE  CRANE  ERROR 


1 .  INTRODUCTION 

1.1  Purpose.  A  concept  named  JUGFAE,  JUG-contained  Fuel-Air  Explosive, 
was  disclosed  (Ref.  1),  tried  (Ref.  2)  and  studied  (Ref.  3)  as  a  method  for 
using  fuel-air  explosions  to  clear  minefields  in  uncontested,  rear  areas. 

The  latter  efforts  have  shown  that  the  concept  has  no  fundamental  flaw. 

Whether  or  not  the  concept  is  practical  can  only  be  found  by  field  experience. 

The  analysis  of  failure  needed  mathematics  that  threatened  to  overrun  the 
concept  report  (Ref.  3)  and  details  were  held  back.  This  report's  Appendix  A 
supplies  the  details  of  the  published  analysis.  Other  appendices  use  separate, 
albeit  agreeable,  methods  that  were  known  to  (Appendix  B)  or  found  afterwards 
by  Umholtz  (Appendix  C,D). 

The  shortened  work  in  Reference  3  dealt  with  the  error  involved  in 
remote  placement  of  numerous  fueled  jugs  on  the  mined  land.  In  Reference  2, 
the  decision  to  use  pre-emplaced  jugs  eliminated  the  remoting  distraction 
from  the  crucial  test  of  detonation  transfer.  That  effort  demonstrated 
that  detonation  wave  transfer  across  overlapping  clouds  would  take  place. 

Previous  methods  (Ref.  3)  arduously  solved  the  failure  probability  for 
a  few  lines  with  up  to  ten  jugs  placed  on  a  line.  The  present  results  give 
the  probability  of  failure  with  crane  placement  for  all  cases  up  to  a  30- jug 
line.  That  length  should  include  all  lines  which  are  practical  to  set  up. 

The  reader  may  only  want  the  results  and  to  know  that  the  six  different 
methods  produce  the  same  low  probability  of  failure  for  the  JUGFAE  concept. 

The  mathematically  inclined  reader  can  profit  by  contrasting  the  methods. 

For  instance,  the  method  described  in  Appendix  D  is  an  achievement  in  concise¬ 
ness  and  new  attack  direction  that  can  only  be  appreciated  after  reading  the 
earlier  appendices.  The  mathematical  efforts  may  be  widely  useful  since  the 
probability  methods  apply  to  a  variety  of  problems. 

1.2  Background .  The  JUGFAE  process  described  crane-emplaced,  fueled 
jugs  spaced  over  mined  land;  one  row  of  the  process  is  depicted  in  Figure  1. 

If  the  carefully  placed  jugs  were  simultaneously  burst,  the  numerous, 
growing  fuel-air  clouds  would  overlap  and,  by  intent,  one  large  explosion 
would  be  triggered  an  instant  later.  Susceptible  land  mines  under  the 
explosion  would  be  neutralized.  The  "don't  cross"  line  would  be  moved  and 
the  setup  repeated.  Advantages  of  the  process  are  greater  safety,  complete¬ 
ness,  and  time  saving  over  one-at-a-time  shooting. 

Overlap  of  the  clouds  is  crucial  to  the  JUGFAE  process  because  a  deton¬ 
ation  wave  cannot  jump  an  air  gap  and  resume  detonating  the  rest  of  the 
clouds.  An  overlap  region  has  to  "chain  together"  every  cloud  to  its  neigh¬ 
bors  in  a  line  or  the  detonation  wave  dies  out.  When  that  happens,  the 
resources  utilized  in  setting  up  the  line  are  partly  wasted  ("partial  failure"). 

Reference  3  cites  two  reasons  why  clouds  might  not  overlap:  a  "critically 
short  cloud"  and  "crane  error".  The  first  reason  was  fully  set  out  there  and 
will  just  be  briefly  reviewed.  The  primary  purpose  of  this  report  is  to  fully 
describe  the  mathematics  of  the  crane  error. 
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2.  PARTIAL  FAILURE 


2.1  Critically  Short  Cloud.  The  first  cause  of  partial  failure  is 
the  inherent  variability  of  the  cloud  growth  process.  Assume  the  jugs  are 
perfectly  placed  on  the  correct  separation  marks.  At  tne  time  when  the  whole 
cloud  row  is  to  be  exploded  (100  ms  after  the  jugs  burst),  some  clouds  will 
overlap  more  and  some  less  than  is  intended.  However,  only  one  disorder 
will  cause  partial  failure.  This  occurs  when  one  cloud  out  of  the  entire  row 
just  touches  its  two  neighbor  clouds.  The  cloud  with  the  critically  short 

radius  (R  . . )  causes  a  tangent  geometry.  The  chaining  together  of  the  clouds 
(represented  by  circles)  is  lost.  The  situation  is  shown  in  Figure  2. 


The  explosion  of  the  entire  cloud  row  is  causec  by  a  single  explosive 
charge,  not  one  charge  per  cloud,  for  reasons  of  simplicity  and  time  saving. 
The  detonation  wave  started  by  the  single  charge  cannot  enter  the  tangent 
cloud  and  so  the  whole  row  is  unexploded.  A  partial  failure  is  recorded 
and  the  cost  is  partly-wasted  time  and  effort  that  went  into  that  setup. 

The  jug  spacing  (S)  is  seen  in  Figure  2  to  be  twice  the  cloud  average 
radius,  shortened  by  a  distance  that  will  statistically  assure  cloud  overlap 
(S  =  2TT  -  PR").  The  amount  of  cloud  overlap  (PIT)  needed  depends  directly  on 
the  variability  of  the  clouds'  radii  (s)  and  on  a  probability  coefficient  (K) 
that  represents  with  stated  confidence  what  portion  of  all  clouds  will  have 
radii  that  will  exceed  the  critically  short  cloud's  radius  (PIT  =  Ks). 

Reference  3  found  the  jug  spacing  for  different  probabilities  that  a  high 
proportion  of  clouds  would  overlap.  For  instance,  clouds  with  a  6.23  m 
average  radius  and  0.25  m  standard  deviation  need  a  jug  separation  of  11.3  m 
(19%  overlap)  to  yield  a  probability  of  99%  that  95%  of  the  clouds  will 
overlap.  (The  probabilities  do  not  depend  on  the  number  of  jugs  in  the  line.) 


DESIRED 


■PR  + 


UNDESIRED 


'TxV 

Rc  r i  t\  /  v 


TANGENT  CLOUD 


Figure  2.  Critically  Short  Cloud. 


2.2  Crane  Placement  of  Jugs.  Figure  3  shows  four  clouds  from  jugs  that 
are  perfectly  placed  and  as  they  could  look  if  the  jugs  were  emplaced  by  a 
crane.  All  the  clouds  are  drawn  as  circles  with  radius  TT.  The  crookedness 
and  mixed  amount  of  overlap  in  the  lower  diagram  are  due  to  crane  error.  The 
crane  sets  a  jug  down  within  an  error  circle  centered  on  the  designated  mark. 
(The  error  circle  bounds  the  miss  distance  from  the  mark  that  we  are  willing 
to  tolerate  or  must  expect  from  crane  emplacement.)  For  statistical  reasons, 
the  diameter  of  the  crane  error  circle  must  be  no  more  than  the  amount  of 
cloud  overlap.  They  are  equal  size  in  Figure  3  and,  as  it  happens,  a  23% 
cloud  overlap  has  been  drawn. 


PERFECT 

PLACEMENT 


CRANE 

ERROR 


Figure  3.  Crane  Placement  Error. 
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One  part  of  Figure  3  also  shows  how  the  tangent  geometry  situation  can 
arise  from  crane  error.  This  very  bad  placement  is  enlarged  in  the  lower 
part  of  Figure  4*  The  point  illustrated  is  that  the  undesirable  tangent 
geometry  is  now  due  to  crane  error,  not  the  critically  short  (radius)  cloud 
described  in  Section  2.1. 

Let  the  diameter  of  the  crane  error  circle  be  2R  .  Let  s  be 

crane  crane 

the  sample  standard  deviation  of  the  miss  distance  of  the  jug  from  the  mark. 
Using  the  sample  standard  deviation,  scrane,  as  an  estimator  of  the  population 

standard  deviation,  <7crane,  set  Rcrane  =  2<rcrane,  i.e.,  the  maximum  tolerable 

miss  distance  of  setting  up  a  jug  is  two  standard  deviations.  For  circular 
symmetry,  and  assuming  the  error  around  the  mark  is  normally  distributed 
(a  uniform  distribution  is  an  alternative  assumption),  it  is  shown  in  Refer¬ 
ence  4  that  the  probability  of  the  event  of  interest,  namely  that  the  jug  is 
outside  the  circle,  is  p  =  0.135.  Otherwise,  the  probability  of  one  jug 
being  set  inside  the  circle  is  q  =  0.Q6 5. 


CRANE  ERROR  CIRCLE  DIA.  =  2R 

CRANE 

JUGS  ON  S=  GOOD  CLOUD  OVERLAP. 
PROBABILITY  ONE  JUG  INSIDE  2  RcranE  =  0  86 


JUGS  ON  OPPOSITE  SIDES 
OF  CRANE  ERROR  CIRCLE, 
TANGENT  CLOUDS. 


IF  JUGS  LAID  ON  S', 
CLOUDS  OVERLAP. 


Figure  4.  Worst  Crane  Placement. 
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An  abstraction  of  Figure  4  shows  that  partial  failure  due  to  crane  error 
proceeds  only  if  several  mistakes  take  place.  First,  the  crane  sets  some 
jugs  outside  the  permitted  error  circle  and  then  these  particular  events 
also  occur:  the  misset  jugs  are  adjacent  ^nd  on  the  left  and  right  sides, 
respectively,  of  the  crane  error  circles.  There  is  no  partial  failure  if 
the  jugs  are  set  any  of  tne  other  three  possible  ways  that  they  can  be  both 
adjacent  and  mislaid.  All  the  possibilities  are  shown  in  Figure  5. 

Note  that  tne  probability  of  failure,  given  tnat  two  adjacent  jugs  are 
mislaid,  is  1/4.  These  four  combinations  are  basic  to  the  analysis  of  crane 
error.  The  idea  illustrated  by  Figure  5  is  used  throughout  this  report. 


CRANE 


TWO  ADJACENT  JUGS  ON 
CRANE  ERROR  CIRCLE. 

LEFT  LEFT 


LEFT  RIGHT*  PARTIAL  FAILURE 


RIGHT  LEFT 


RIGHT  RIGHT 


Figure  5.  Misplacements  of  Adjacent  Jugs. 


Left  and  right  are  just  convenient  names  if  one  is  next  to  the  crane, 
watching  it  set  in  jugs.  Left  and  right  do  not  refer  at  all  to  which  way 
the  detonation  wave  is  moving. 
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3.  SUMMARY  ANALYSES  OF  CRANE  ERROR 


3.1  Diagrams .  The  original,  unpublished  solution  to  the  crane  error 
problem  is  due  to  J.  Richard  Moore,  Ballistic  Research  Laboratory,  retired. 

It  depended  on  conditional  probability  and  seldom  used  mathematics  of  the 
probability  of  runs.  The  arcane  run  portion  of  the  solution  was  replaced 

in  Reference  3  with  the  diagrammatic  method  of  counting  the  space  of  events. 
Although  results  for  longer  lines  were  given,  the  only  worked-out  example 
was  for  a  2- jug  line.  The  diagrammatic  method  is  shown  for  several  cases 
in  Appendix  A;  Moore's  run  solution  is  in  Appendix  B. 

The  basic  equation  used  to  calculate  failure  probability  is  given  in 
Appendices  A  and  B.  It  is 

P(failure)  =  ^^P(fji  out)P(i  out).  (1) 

i=l 

To  calculate  P(failure;  using  Equation  1,  take  the  probability  of  one 
jug  of  the  line  of  n  being  out  [P(l  out)]  and  multiply  by  the  probability  of 
failure  for  that  arrangement,  i.e.,  of  one  jug  being  out  of  the  crane  error 
circle  [P(f  J 1 ) ] ;  taka  the  probability  of  two  jugs  of  the  line  of  n  being  out 
[P(2  out)]  and  multiply  by  the  probability  of  failure  for  that  arrangement, 
i.e.,  of  two  jugs  being  out  of  the  crane  error  circle  [P(f|2)].  Continue  to 
form  the  product  of  these  two  terms  all  the  way  through  the  case  of  n  jugs 
being  mislaid.  Add  all  the  products  and  the  sum  is  the  probability  of 
failure  -for  an  n-jug  line. 

The  method  presented  in  Appendix  A  is  to  make  individual  diagrams  of 
one  jug  mislaid  of  the  n-jug  line,  two  jugs  mislaid  of  the  n-jug  line,  ...  , 
n  jugs  mislaid  of  the  n-jug  line.  Since  the  line  cannot  fail  with  only  one  jug 
out,  i.e.,  P ( f | 1  out)  =  0,  no  diagram  needs  to  be  drawn.  For  two  jugs  out, 
all  the  diagrams  must  be  drawn.  A  diagram  with  two  adjacent  jugs  mislaid 
must  be  broken  down  further  because  only  a  left-right  misplacement  on  the 
crane  error  circle  leads  to  detonation  failure  (Section  2.2).  The  diagrams 
are  counted  for  failure  combinations  and  the  count  quickly  gives  the  term 
P ( f ! i  out).  Since  the  number  of  diagrams  grows  exponentially,  e.g.,  a  6-jug 
line  needs  2**6  or  64  diagrams,  another  technique  is  needed  to  find  P(f|i>6). 

3.2  Runs  of  Jugs.  The  method  presented  in  Appendix  B  starts  with  Equa¬ 
tion  1  but  differs  from  Appendix  A  by  finding  P(f|i  out)  utilizing  run  theory 
instead  of  diagrams.  (The  term  is  calculated  by  Equation  B2  in  a  fashion 
similar  to  Equation  1).  Equation  B2  treats  i  jugs  being  out  of  place  as 
meaning  runs  of  misplaced  jugs  have  formed.  A  run  is  an  unbroken  sequence. 

One  mislaid  jug  is  a  run  of  one.  Two  mislaid  jugs  could  mean  two  adjacent 
jugs  are  mislaid  or  two  jugs  anywhere  on  the  line  are  mislaid;  the  former 

is  one  run  of  length  two  and  the  latter  is  two  runs  of  length  one.  Each  run 
that  can  be  formed  for  i  jugs  being  mislaid  has  a  probability  of  occurrence 
that  can  be  calculated  using  run  theory,  (see  Equation  B3).  Each  run  also 
has  a  probability  of  failure,  e.g.,  probability  of  failure  of  one  run  of 
length  two  is  1/4  as  shown  in  Figure  5*  The  two  probabilities  are  multiplied; 
all  runs  with  mislaid  jugs  amounting  to  i  are  so  treated  and  the  products 
are  summed  to  give  P(f|i  out).  The  term  P(i  out)  is  found  by  the  binomial 
distribution  (Equation  A3)  from  which  Equation  1  is  then  computed.  The  theory 
of  runs  is  clarified  in  Appendix  B.  Approximations  are  introduced  to  save 
labor,  however,  and  cause  a  misleading  indication  that  the  probability  of 
failure  levels  off  after  six  or  ten  jugs;  it  actually  rises  almost  linearly 
with  jug  number.  We  have  not  gone  past  30  jugs,  but  obviously,  for  some 
large  n,  failure  of  the  line  is  almost  certain.  Exact  calculations  are 
reported  in  Appendices  C  and  D. 
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3.3  Runs  of  Arrangements.  Moore's  run  method  combined  the  conditional 
probability  of  a  jug  being  misplaced  with  the  probability  of  a  run  of  such 
jugs.  Conditional  probability  is  not  used  in  Appendix  C  but  counting  of 
runs  is.  The  general  idea  is  as  follows.  All  possible  runs  of  misplaced 
jugs  are  listed  for  an  n-jug  line.  The  product  of  three  facLors  determines 
the  probability  of  failure  of  a  run.  All  the  run  probabilities  are  summed 
to  yield  the  probability  of  failure  with  an  n-jug  line.  Examples  for  3-  and 
5- jug  lines  are  given  in  Appendix  C. 

The  idea  in  more  detail  is  as  follows.  An  n-jug  line  is  regarded  in 
Appendix  C  as  an  arrangement  of  jugs,  whether  outside  or  inside  of  the  error 
circle.  If  all  jugs  are  correctly  placed,  this  arrangement  is  described  as 
one  run  of  length  n.  There  is  a  probability  of  that  arrangement  happening 
and  there  is  only  one  such  run  possible.  The  maximum  number  of  runs  is  the 
integer  part  of  (n+l)/2  for  odd  lines  and  n/2  for  even  lines,  e.g.,  a  5- jug 
line  has  a  maximum  of  three  runs.  (The  maximum  number  of  runs  occurs  in  the 
jug  arrangement  where  one  end  jug  is  out,  the  next  jug  is  in,  the  next  out, 

...»  and  tne  other  end  jug  is  out.)  Each  run  can  occur  in  different 
arrangements  along  the  line,  e.g.,  a  run  of  one  in  a  5- jug  line  can  occur 
in  five  ways.  In  Appendix  C,  all  runs  from  one  to  the  maximum  and  all 
arrangements  having  that  number  of  runs  are  found.  The  number  of  arrangements 
can  be  diagrammed  or  calculated  by  Equation  C2.  The  probability  of  getting 
so  many  jugs  out-and-in.  Equation  Cl,  the  probability  of  failure  of  a 
particular  run,  Equations  C4  and  C5»  and  the  number  of  arrangements  possible, 
Equation  C2,  are  multiplied  together  and  the  individual  run  probabilities  are 
summed  to  get  the  probability  of  failure  with  an  n-jug  line,  Equation  C7. 

The  procedure  is  programmed  in  "JUGRUN". 

3.4  The  Bad  Event.  The  solution  of  the  crane  error  problem  described 
in  Appendix  D  is  distinctly  different  from  the  related  approaches  of  the 
three  previous  appendices.  The  solution  does  not  use  either  run  theory  or 
conditional  probability.  Instead  it  recognizes  the  bad  event  of  Figure  4  as 
one  that  can  arise  n-1  times  in  an  n-jug  line.  The  bad  event  E  occurs  when 
adjacent  jugs  are  set  on  the  left  side  and  the  right  side  of  their  respective 
crane  error  circles.  Let  event  E,,  be  so  setting  jug  1  and  jug  2;  event  Ep  is 
so  setting  jug  2  and  jug  3;  ...  ; 'event  Eq_1  is  so  setting  jug  (n-1 )  and  jug  n. 
Failure  (partial  detonation  of  the  n-jug  line)  will  occur  if  at  least  one  of 
the  (n-1)  events  i3  realized. 

P (failure )  =  P(E1  U  E2  U  E3  \j  . . .  U  En_1 ) ,  (2) 

where  "y"  symbolizes  "or". 

The  general  expansion  of  Equation  2  is  given  in  Appendix  D  and  is 
performed  two  ways,  leading  to  two  dissimilar  equations,  D5  and  D25,  for 
the  probability  of  failure.  The  two  ways  are  indicated  simply  by  the 
expansion  of  the  union  of  '  y  two  3ets  A  and  B, 

A  U  B  '.  +  B  -  (A  H  B),  (3) 

A  U  B  =  A  (B  O  A),  (4) 

where  "p"  symbolizes  "and",  "A"  symbolizes  "not  A". 

Venn  diagrams  representing  this  union  are  shown  in  Appendix  D.  An  important 
fact  used  for  simplification  in  both  approaches  is  that  the  intersection  of 
adjacent  events  yields  the  null  set  because  they  cannot  occur.  For  instance, 
in  order  for  E.  u  Ep  to  occur,  it  would  require  jug  2  to  be  both  too  far  right 
and  too  far  left. 
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The  full  expansion  of  Equation  2,  in  the  manner  of  Equation  3,  is 
presented  in  Appendix  D.l.  Most  of  the  terms  of  Equation  D2  will  go  to 
zero  because  of  the  intersection  of  mutually  exclusive  events.  This 
simplification  will  be  illustrated  by  the  example  of  a  5-jug  line.  The 
remaining  terms  for  an  n-jug  line  are  represented  by  Equation  5.  To  avoid 
calculating  a  factorial  term,  a  recursive  formula  is  given  in  Equation  D7 
that  is  better  suited  to  computation.  The  procedure  is  programmed  in 
"JUGCOM". 


P(failure) 


(5) 


where , 
n  = 
In/2J  = 
P  = 
p/2  = 

P'  = 


number  of  jugs  in  a  row, 
the  integer  part  of  n/2, 

probability  of  an  individual  jug  being  outside  error  circle, 

probability  of  an  individual  jug  being  too  far  left  OR 

ility  of  an  individual  jug  being  too  far  right, 

=  probability  of  one  jug  being  too  far  left  and  the  next  jug 
being  too  far  right  (by  independence  of  placements  of  jugs)  OR 
probability  of  a  particular  pair  of  consecutive  jugs  causing  failure. 


probab 

(p/2; 


The  full  expansion  of  Equation  2,  in  the  manner  of  Equation  4»  is  in 
Appendix  D.2.  A  2-  and  a  3- jug  line  are  simply  expressed,  but  a  4- jug  line 
requires  numerous  set  manipulations  before  arriving  at  the  probability  of 
failure  as  the  sum  of  the  failure  probabilities  for  the  3- jug  line  and  a 
constant  times  the  failure  probability  of  the  2- jug  line  (Equation  D18). 

The  identical  form  is  given  for  the  n-jug  line  by  Equation  6.  The  procedure 
is  programmed  in  "JUGREC". 


P(failure)  =  P(fail  w/n-1  jugs)  +  p'[l  -  P(fail  w/n-2  jugs)].  (6) 


In  Appendix  D.3,  an  approximate  formula  for  the  probability  of  failure 
is  given  by  Equation  D30.  The  derivation  starts  with  the  general 
Equation  2,  but  wrongly  assumes  that  the  adjacent  E  events  can  occur. 

The  deliberate  error  is  made  in  Equation  D28.  The  approximation  is  quite 
good  for  any  jug  line  length  and  using  probabilities  of  missetting  one  jug 
around  p  ^  0.135. 

2  n-1 

P(failure)  =*  1  -  [1  -  (p/2)  ]  .  (7) 
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4.  RESULTS 


The  analyses  yielded  the  probability  of  failure  to  detonate  overlapped 
clouds  from  a  row  of  jugs  set  in  by  a  crane.  Table  1  was  calculated  with 
a  probability  of  missetting  one  jug  of  0.135,  for  reasons  set  forth  in 
Section  2.2,  and  the  failure  probability  for  up  to  a  30-jug  row  is  calculated. 
Table  2  presents  the  failure  probability  of  the  row  of  jugs  for  other  proba¬ 
bilities  of  missetting  one  jug  around  the  value  of  0.135.  Figure  6  is  a 
linear  fit  of  the  data  in  Table  2. 
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TABLE  1.  Probability  of  Failure  Due  to  Crane  Error. 


Number 

Probability 

of  Jugs 

of  Failure 

in  Row 

Due  to  Crane 
Error  (*) 

1 

0.0000 

2 

.0046 

3 

.0091 

4 

.0136  . 

5 

.0182 

6 

.0227 

7 

.0271 

8 

.0316 

9 

.036 0 

10 

.0404 

11 

.0448 

12 

.0492 

13 

•  .  .0535 

14 

.0579 

15 

.0622 

16 

.0665 

•17 

.0708 

18 

.0750 

19 

.0792 

20 

.0835 

21 

.0877 

22 

.0918 

23 

.0960 

24 

.1001 

25 

.1042 

26 

.1083 

27 

.1124 

28 

.1165 

29 

.1205 

30 

.1246 

(*)  The  crane  is  permitted  to  set  the  jug  off  the  ideal  mark  but 
within  an  error  circle  whose  radius  is  equivalent  to  half  the  cloud  overlap 
distance,  P K.  The  probability  of  any  individual  jug  being  outside  the  circle 
is  0.135.  The  probability  of  failure  due  to  a  short  radius  cloud  is  about  0.01 
(by  Section  2.1  and  Reference  3).  The  total  probability  of  failure  due  to 
either  a  short  radius  cloud  or  crane  error  is  found  by  adding  0.01  to  the 
numbers  in  column  two. 
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TABLE  2.  JUGFAE  Sensitivity. 

Jugs  Probability  of  Missetting  One  Jug,  (p). 

p=.1  p  =  .135  p=.15  p=.2 

n  Probability  of  Failure  for  n-Jug  Line 


2 

0.0025 

0.0046 

0.0056 

0.0100 

3 

.0050 

.0091 

.0113 

.0200 

4 

.0075 

.0136 

.0168 

.0299 

5 

.0100 

.0182 

.0224 

.0397 

6 

.0125 

.0227 

.0279 

.0494 

7 

.0149 

.0271 

.0334 

.0590 

8 

.0174 

.0316 

.0389 

.0685 

9 

.0199 

.0360 

.0443 

.0779 

10 

.0223 

.0404 

.0497 

.0872 

11 

.0248 

.0448 

.0551 

.0965 

12 

.0272 

.0492 

.0605 

.1056 

13 

.0297 

.0535 

.0658 

.1146 

14 

.0321 

.0579 

.0711 

.1236 

15 

.0345 

.0622 

.0763 

.1324 

16 

.0369 

.0665 

.0815 

.1412 

17 

.0393 

.0708 

.0867 

.1499 

18 

.0418 

.0750 

.0919 

.1584 

19 

.0442 

.0792 

.0970 

.1669 

20 

.0466 

.0835 

•  .1022 

.1754 

21 

.0489 

.0877 

.1072 

.1837 

22 

.0513 

.0918 

.1123 

.1919 

23 

.0537 

.0960 

.1173 

.2001 

24 

.0561 

.1001 

.1223 

.2082 

25 

.0584 

.1042 

.1273 

.2162 

26 

.0608 

.1083 

.1322 

.2241 

27 

.0632 

.1124 

.1371 

.2319 

28 

.0655 

.1165 

.1420 

.2397 

29 

.0678 

.1205 

.1468 

.2474 

30 

.0702 

.1246 

.1517 

.2550 

12 


Probability  of  Failure  x  100 


0  6  10  16  20  26  30 


n-Jug  Line 

Figure  6.  Linear  Plot  of  Failure  to  Detonate  the  Line. 
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5.  CONCLUSION 


Using  the  JUGFAE  concept,  an  uncontested  minefield  could  be  cleared, 
a  line  at  a  time,  by  one  detonation  of  overlapping  fuel-air  clouds.  The 
clouds  are  produced  from  fuel  jugs  set  down  by  a  crane  operating  parallel  to 
and  alongside  the  edge  of  the  minefield.  If  the  jugs  could  be  ideally  spaced, 
the  vagaries  in  cloud  size  would  not  prevent  their  necessary  overlap.  The 
crane  will,  however,  set  each  jug  around  the  ideal  nark.  This  crane  error 
can  prevent  cloud  overlap  unless  each  jug  is  set  down  within  the  allowable 
error  circle  around  the  ideal  mark. 

This  report  presents  different  analyses  of  the  probability  of  failure, 
caused  by  crane  error,  to  detonate  clouds  from  a  line  of  jugs.  The  methods 
used  include  numeration  of  event  space,  conditional  probability,  run  theory 
and  set  theory,  either  separately  or  combined.  Set  theory  proved  to  be  much 
shorter,  both  conceptually  and  computationally,  than  the  other  methods. 

One  result  is  that  given  an  expected  probability  of  0.135  of  missetting 
one  jug,  the  probability  of  not  detonating  an  entire  10-jug  line  is  0.0404. 
Furthermore,  the  probability  of  failure  is  a  little  less  than  proportional  to 
the  number  of  jugs  in  the  line.  The  major  result  is  that  crane  error  results 
in  a  manageably  small  probability  of  failure  to  detonate  the  entire  line. 

The  extensive  analysis  supports  the  conclusion  that  there  is  no  basic  flaw 
with  the  JUGFAE  concept.  However,  the  concept's  operational  feasibility  has 
yet  to  be  tested. 
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APPENDIX  A 

DIAGRAMS  OF  CRANE  ERROR 
FOR  THREE-,  FOUR-,  AND  FIVE-JUG  LINES 
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A.1  Explanation,  The  harm  in  setting  the  jugs  off  the  spacing  mark  is 
the  likely  nonoverlap  of  some  clouds  and  subsequent  failure  to  propagate 
the  detonation.  The  probability  of  failure  of  a  two-jug  line  was  analyzed 
completely  in  Reference  3.  Diagrams  showing  all  possible  placements  were 
used.  The  method  is  extended  here  to  more  lines.  Definitions  of  events 
leading  to  failure  and  the  probability  of  failure  equation  are  reviewed. 

The  calculations  and  respective  diagrams  are  listed  for  three-,  four-,  and 
five-jug  lines.  The  diagrams  use  (+)  for  a  jug  on  the  mark  and  (0)  for  a  jug 
set  outside  the  permitted  error  circle  around  the  mark.  A  failure  situation 
(*)  occurs  when  adjacent  jugs  are  placed  on  the  left  (L)  and  right  (R)  sides 
of  their  respective  error  circles,  producing  tangent  clouds.  (See  Figure  4.) 

For  a  cloud  line  detonation  failure  to  occur,  three  events  must  take  place 
during  jug  placement: 

Eq  =  the  crane  puts  2  or  more  jugs  outside  their  permitted  error  circles, 

E^  =  given  that  Eg  has  occurred,  at  least  two  outside  jugs  are  adjacent, 

=  given  that  E.  has  occurred,  adjacent  jugs  are  on  the  left  and  right 
sides,,  respectively,  of  their  error  circles. 

Failure  to  transmit  detonation  with  crane  placement  of  jugs  can  be  restated 
as:  failure  occurs  if,  given  that  EQ  has  occurred,  then  E1  and  E2  occur. 

The  probability  of  failure  P(failure)  can  be  written 

n 

P(failure)  =  L  Pi(E0)  Pi(E2,E1|E0).  (A1) 

i=0 

The  first  term  is  found  by  the  binomial  distribution, 

Pi(E0)  =  cr  P^1  ~  P)n~r»  r  =  (A2) 

n  =  total  number  of  jugs 

On! 

=  -  ,  n!  =  n(n-l)(n-2)...(3)(2)(l).  (A3) 

r!  (n-r) ! 

The  second  term,  a  conditional  probability  term,  is  found  using  diagrams 
which  will  be  listed  later  in  this  Appendix.  Equation  A1  will  be  expanded 
for  the  three  cases  to  be  discussed. 
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Three- jug  line. 

P(failure)  =  P(0  out)P(E2>E1  jO  out)  +  P(l  out)p(E2, E-,  { 1  out) 

+  P(2  out)P(E2,E1 j2  out)  +  P(3  out)P(E2, E1 j 3  out). 


(A4) 


From  Section  2.2,  the  simple  probability  of  mislaying  one  jug,  i.e.,  setting 
it  outside  the  crane  error  circle,  is  p  =  0.135.  By  Equations  A2  and  A3, 

3  0  3-0  0  3 

P  (E  )  =  P(0  out)  =  C  p  (1-p)  =  1 ( .135)  (.865)  =  0.6472 

0  0  0 

313-1  12 

P  (E  )  =  P(1  out)  =  C  p  (1-p)  =  3( . 1 35)  (.865)  =  0.3030 

1  0  1 

323-2  21  (A5) 

P  (E  )  =  P(2  out)  =  C  p  (1-p)  =  3 ( . 1 35 )  (.865)  =  0.0473 

2  0  2 

333-3  30 

P  (E  )  =  P(3  out)  =  C  p  (1-p)  =  1 ( . 135)  (.865)  =  0.00246 

3  0  3 

The  conditional  probabilities  are  brought  forward  from  the  3- jug  diagram. 
P(failure)  =  0.6472(0)  +  0.3030(0)  +  0.0473(1/6)  +  0.00246(4/8), 

P(failure)  =  0.00911  for  a  3- jug  line  as  shown  in  Table  1.  (A6) 


Four- jug  line. 


P(failure)  = 

=  P(0  out)P(E2, 

Eii 

0  out) 

+  •  •  •  H 

1-  P(4  out)P(E2,E1  i4  out). 

(A7) 

P(0 

OUt)  : 

=  C4p°(l-p)4  = 

K 

.135)°( 

.865)4 

=  0.560 

P(1 

OUt)  : 

=  c^d-p)3  = 

4( 

-135)1( 

.865) 3 

=  0.349 

P(2 

OUt)  : 

=  c4p2(l-p)2  = 

6( 

- 1 35 ) 2  ( 

.865)2 

=  0.082 

(A8) 

P(3 

OUt)  = 

=  C^d-p)1  = 

4( 

•135)3( 

.865)1 

=  0.0085 

P(4 

out)  : 

II 

O 

4^  -fx 

X) 

4*. 

—X 

1 

0 

1( 

•135)4( 

.865)° 

=  0.00033 

The  conditional  probabilities  are  brought  forward  from  the  4- jug  diagram. 
P(fail)  =  0.560(0)  +  0.349(0)  +  0.082(1/8)  +  0.0085(12/32)  +  0.00033(11/16), 
P(fail)  =  0.0136  for  a  4-jug  line  as  shown  in  Table  1.  (A9) 

0 
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Five- jug  line. 


P(failure ) 
P(0  out) 

P(1  out) 

P(2  out) 

P(3  out) 

P(4  out) 

P(5  out) 


P(0  out)P(E2fE1 |o  out)  +  ...  +  P(5  out)P(E2,E1 

5  0  5  0  5 

C  p  (1-p)  =  1(.135)  (.865)  =  0.4842 
0 

5  14  14 

C  p  (1-p)  =  5(.135)  (.865)  =  0.3779 


1 

5  2  3 

C  p  (1-p) 

2 

5  3  2 

C  p  (1-p) 

3 

5  4  1 

C  p  (1-p) 

4 

5  5  0 

C  p  (1-p) 

5 


2  3 

10( .135)  (.865)  =  0.1179 


3  2 

1 0 ( .135)  (.865)  =  0.0184 


4  1 

5( .135)  (.865)  =  0.00144 


5  0 

1 ( .135)  (.865)  =  0.00004 


5  out) . 


(A10) 


(All) 


The  conditional  probabilities  are  brought  forward  from  the  5- jug  diagram. 

P(failure)  =  0.4842(0)  +  0.3779(0)  +  0.1179(1/10)  +  0.0184(24/80) 

+  0.00144(45/80)  +  0.00004(26/32). 

P(failure)  =  0.01815  for  a  5-jug  line  as  shown  in  Table  1.  (A12) 

A. 2  Diagrams.  All  the  arrangements  of  three-,  four-,  and  five-jug  lines 
are  drawn  in  the  way  explained  in  Section  A.1.  The  equations  of  Section  A.1 
help  account  for  all  diagrams.  The  number  of  arrangements  is  Cn  and 
since  each  jug  out  (up  to  r)  can  be  mislaid  left  and  right,  there  are  2r 
sub-pictures  under  each  arrangement.  Not  all#diagrams  need  to  be  drawn 
explicitly.  The  mirror  images  are  redundant.  (With  an  arrangement  read 
from  left  to  right,  its  mirror  image  is  the  one  read  right  to  left.  For 
example,  +000  and  000+  are  mirror  images;  the  probability  of  failure 
is  4/8  for  each. ) 


A  reason  to  draw  all  diagrams  is  that  they  directly  inform  a  project 
officer  of  the  probability  of  failure  if  he  has  knowledge  of  how  many  or 
which  particular  jug3  are  mislaid. 


Dixon,  Lisa  A.,  personal  communication,  Aberdeen  Research  Center,  Aberdeen, 
Maryland,  August  1988. 
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THREE-JUG  LINE 


no  jugs  out 


3 

C  =  1 
0 

0 

2  =  1 


3 

C  =  3 
1 


2=2 


3 

C  =  3 
2 

2 

2=4 


3 

C  =  1 
3 

3 

2=8 


+ 

+ 

+ 

0/0 

one 

jug 

out 

+ 

+ 

_ 0 

+ 

+ 

L 

+ 

+ 

R 

0/0 

p(e. 

> » E^  ! 

[1  out) 

two 

jugs  out 

+ 

0 

_ 0 

+ 

L 

L 

+ 

L 

R  * 

+ 

R 

L 

+ 

R 

R 

1/4 

p(e 

2,E1 

|2  out) 

three  jugs  out 

0 _ 

0 

0 

L 

L 

L 

L 

L 

R  * 

L 

R 

L  * 

L 

R 

R  * 

R 

L 

L 

R 

L 

R  * 

R 

R 

L 

R 

R 

R 

P(E2’®1 ' ^  =  0 


0/0 

(0  +  0  +  0)/3  =  0 


=  (1/4  +  0  +  l/4)/3  =  1/6 


4/3 

P(E2,E1|3  out)  =  4/S 
*  failure  situation 
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FOUR-JUG  LINE 
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P(E2,E1 ! 3  out)  =  (4/3  +  2/8  +  2/8  +  4/8 )/4  =  12/32 


*  failure  situation 


23 


FOUR-JUG  LINE  (continued) 
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P(E  e1 1 4  out)  =  11/16 
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FIVE-JUG  LINE  (continued) 
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FIVE-JUG  LINE  (continued) 
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*  failure  situation 
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FIVE-JUG  LINE  (continued) 


four  jugs  out 


LLL  +  L  LLLL  + 

LLL  +  R  LLLR+* 

llr+l*llrl+* 

L  L  R  +  R  *  L  L  R  R  +  * 

lrl+l*lrll+* 

lrl+r*lrlr+* 

lrr+l*lrrl+* 

lrr+r*lrrr+* 

RLL  +  L  RLLL  + 

RLL  +  R  RLLR+* 

R  L  R  +  L  *  R  L  R  L  +* 

RLR*R*RLRR+* 

RRL  +  L  RRLL  + 

RRL  +  R  RRLR+* 

RRR  +  L  RRRL  + 

RRR  +  R  RRRR  + 

8/16  11/16 

P(E2,E1i4  out)  =  (11/16  +  8/16  +  7/16  +  8/16  +  11/l6)/5  =  45/60 


*  failure  situation 
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FIVE-JUG  LINE  (continued) 
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26/32 

P(E2,E1 J 5  out)  =  26/32 


*  failure  situation 


five  jugs  out 
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L  L  L  L  L 

L  L  L  L  R  * 

L  L  L  R  L  * 

L  L  L  R  R  * 
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APPENDIX  B 

THEORY  OF  RUNS  APPLIED  TO  CRANE  ERROR 
FOR  THREE-  AND  TEN-JUG  LINES 
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(B1) 


The  probability  of  failure  due  to  jug  placement  error  by  the  crane  is 

n 

P(failure)  =  ^  P  .  (F.n)P .  (E?  ,F. ,  j  EQ)  , 
i=0 

where  the  E  events  were  defined  in  Appendix  A.  The  first  term  is  the  proba¬ 
bility  of  the  crane  placing  i  jugs  outside  the  error  circles.  The  term  is 
calculated  from  the  binomial  distribution  (see  Equation  A3).  The  second 
term  is  the  probability  of  adjacent  jugs  being  on  the  left  and  right  sides  of 
their  respective  error  circles  (Figure  4)*  The  theory  of  runs  can  replace 
the  diagrams,  to  a  large  extent,  for  calculating  the  second  term. 

pi(E2»E1|Eo)  =  P(E2,E1  Srun)P(run'i  outside  error  circle).  (B2) 

all  possible  runs 

The  compound  event  E^E^  is  failure  to  detonate  the  whole  line  of  clouds,  i.e., 
"partial  failure". 

A  run  is  an  unbroken  sequence,  which  in  this  application  refers  to  a  jug 
placed  outside  the  error  circle  of  the  crane  (0).  A  jug  being  on  the  ideal 
mark  (+). starts  a  run  of  another  kind.  The  runs  for  the  +'s  do  not  concern  us. 
Consider  a  3-jug  line  and  the  sequences  that  can  occur  with  one  jug  out. 

+  +  0  +0  +  0  +  + 


There  are  three  runs  of  length  one  for  the  0's  and  they  are  designated  1^, 
i.e.,  (run  lengthnumber  Qf  runs),  designated  by  rfc.  (B3) 

With  two  jugs  out,  the  sequences  are  drawn  as: 

+  00  0  +  0  00  + 


Obviously  the  probability  of  a  run  of  length  two  is  2/3.  For  even  low  numbers 
of  jugs  (see  Appendix  A),  it  is  difficult  to  picture  all  the  ways  the  jugs 
could  be  misplaced.  The  theory  calculates  the  probability  of  a  particular 
run  arising.  From  Reference  B1 ,  the  probability  of  having  exactly  k  runs  of 
0's  amounting  to  r^  in  all  is, 


P(rk) 


C:  :)C;  VC?) 


and  the  probability  of  the  0's  being  arranged  in  v  multiple  runs  is, 


rr 


(B4) 


(B5) 


where  r^  is  the  number  of  0's  and  r2  is  the  number  of  +'s;  k^+  k2+...+  ky=  k. 


Feller,  William,  An  Introduction  to  Probability  Theory  and  Its  Applications, 
Vol  I,  Third  Edition,  John  Wiley  &  Sons,  Inc.,  p.  62. 
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The  notation  used  in  Equations  B4  and  B5  is  the  familiar. 


(I) 


n! 


=  -  ,  where  r  n  are  any  positive 

r!(n-r)'  integers. 


(B6) 


For  the  3- jug  line  with  two  jugs  out  of  place  there  are,  from  the  drawing 
after  Equation  B3,  only  two  possible  sequences:  one  run  of  length  two  or  two 
runs  of  length  one.  So  for  run  2^;  r.,  =  2,  k  =  1 ;  r9  =  3  -  r.,  =  1 


1 


1 


r.  =2,  k=2;r9=3-r.  =1 


From  Equations  B3  and  B4  the  probability  of  the  single  run  arising  is, 

P(21 

(For  a  3-jug  line  the  run  2^  looks  like:  .+  0  0  or  0  0  +  ) 


■•■-ex  D/rn-oo/o- 


-v-c:x;j/r)-- 


(B7) 


(For  a  3- jug  line  the  run  12  looks  like:  0  +  0  ) 

Equation  B5  also  gives  P(l2)  when  k  =  2,  k1  =  2,  k2  =  k?  =  0. 

Note  that  the  run  probabilities  agree  with  the  probabilities  from  the 
explicit  pictures  in  Appendix  A  and  also  that  the  probabilities  sum  to  unity. 


The  first  term  in  Equation  B1  can  be  determined  by  the  condition  for 
failure,  i.e.,  a  left-right  jug  placement  occurs,  by  using  common  sense  or  the 
pictures  in  Appendix  A.  For  example,  the  run  12  cannot  cause  failure  since 
there  are  not  two  adjacent  jugs.  If  the  run  2^  occurs  (the  probability  is 
2/3  that  it  will,  if  the  crane  mislays  two  jugs  out  of  three)  failure  is 
possible.  The  probability  of  failure  is  l/4  when  the  run  does  occur  because 
of  the  simple  arrangement  +00  or  for  the  00+  arrangement. 

+  L  L 

+  L  R  *  (failure) 

+  R  L 

+  R  R  1/4 

The  3-jug  example  will  now  be  written  out.  In  place  of  the  event  notation, 
we  write:  f  =  failure  =  E^,e2,  so  Equation  B1  restated  for  a  3-jug  line  is: 
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Three- jug  line: 

P(failure)  =  P(0  out)P(f|0  out)  +  P(1  out)P(fjl  out) 

+  P(2  out)P(f|2  out)  +  P(3  out)P(f|3  out). 

The  conditional  probabilities  require  Equation  B2  to  be  restated, 

P(f|0  out)  =  0 

P(f|l  out)  =  P(l1)p(fjl1)  =  (1)(0)  =  0 

P(f|2  out)  =  P(l2)p(f‘l2)  +  P(21)P(f|21)  =  (l/3)(0)  +  (2/3)(l/4)  =  2/12 

P(f i 3  out)  =  P(31)p(f j31)  =  (l)(l/2)  =  1/2 

Recall  from  Appendix  A,  the  3- jug  line,  Equation  A6, 

P(0  out)  =  C3p°(i  _  p)3-0  =  C^(o.135)°(0.865)3  =  0.6472,  etc. 

Now  Equation  B8  can  be  evaluated, 

P(failure)  =  (0.6472)(0)  +  0.3030(0)  +  0.0473(2/12)  +  0.00246(1/2) 

P(failure)  =  0.00911  for  three  jugs,  in  agreement  with  Appendix  A  and  Table  1. 

Ten- jug  line; 

Equation  B1  with  failure  (f)  =  (E^.E^)  for  ten  jugs  is  written: 

P(f)  =  P(0  out)P(f 1 0  out)  +  P(1  out)P(f *  1  out)  +  ...  +  P(10  out)p(f|l0  out). 

Since  events  (E^jE-j)  cannot  occur  with  zero  or  one  jug  mislaid, 

P(f jO  out)  =  P(fjl  out)  =0.  A  further  simplification  by  the  binomial 
distribution  (Equation  A3)  is  that  P(i  out)  <  10E-3  for  i  >  6,  i.e.,  there 
is  not  much  chance  that  the  crane  will  mislay  six  jugs.  The  calculation  of 
Equation  B9  contains  a  product  involving  p  =  (0.135)E6  o:  10E-5>  so  the 
equation  can  be  simplified  to  a  much  shorter  summation. 

P(f)  =  P(2)P(f 1 2  out)  +  P(3)P(f j  3  out)  +  P(4)P(f'4  out)  +  P(5)P(f|5  out)  +  e, 
where  e  <  10E-3. 

The  conditional  probabilities  are,  by  Equation  B2,  all  combinations  of  runs 
which  would  place  i  jugs  outside  the  crane  error  circle. 

P(fj2  out)  =  P(l2)P(fjl2)  +  P(21)P(f|21)  =  (36/45)(0)  +  (9/45)0/4)  =  0.0500, 
by  Equation  B4  and  the  simple  diagram  of  probability  of  failure  for  two  jugs. 

P(f|3  out)  =  P(l3)P(f|l3)  +  P(l1,21)P(fjl1,21)  +  P(31)P(f{31) 

=  ( 56/120)  (0)  +  (56/120)  (1/4)  +  (8/120)0/2)  =  0.150. 


(B8) 


(B9) 


(B10) 


(Bll) 
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Note  that  a  multiple 


P(1 


1 ,21 )  - 


run  occurred  in  Equation  B11  so  Equation  B5  is  needed. 


(B12) 


The  probability  of  failure  for  the  multiple  run  (1^,2*)  is  argued  as 
follows.  For  success,  the  detonation  wave  must  travel  through  one  of  the 
runs,  say  '\^t  then  travel  through  tne  next  run  2^.  Since  the  two  events  are 
independent  of  each  other,  the  probability  of  doing  both  is  the  product 
B(l-|  )P(2^ )  {success;  and  tne  probability  of  not  traveling  the  10-jug  line 

holding  these  particular  runs  is  the  complement. 


P(fil1,21)  =  1  -  P(11)P(21) isuccess  =  1  -  (1 ) (3/4)  =  1/4- 
P(f,'4  out)  =  P(14)P(f'i4)  +  P(41)P(f|41)  +  (B13) 

+  P(22)P(f|22)  +  P(12,21)P(f',12,21) 


=  (35/210)(0)  +  (7/210) (ll/l6)  +  (42/210) [l  -  (1 ) (4/8)3 
+  (21/210) [1  -  (3/4) **2]  +  (105/210)(l/4) 

=  0.2917. 

P(fi5  out)  =  P(l5)p(f|l5)  +  P(l3,21)P(fH3,21)  +  P(l2,31)P(fil2.3i) 

+  P(11,22)P(f|l1,22)  +  P(21,31)P(f|2v31) 

+  P(11,41)P(fH1,41)  +  P(51)P(f|51) 

=  (6/252)(0)  +  (60/252) (l/4)  +  (60/252) (l/2) 

+  (60/252) [1  -  (3/4) **2]  +  (30/252) [l  -  1(11/16)] 

+  (30/252) [1  -  (3/4)(l/2) ]  +  (6/252)(26/32) 

=  0.4583. 


Note  that  the  runs  4-j  and  5^  use  the  diagrams  of  Appendix  A  to  supply  the 
conditional  probabilities  P ( f ' 4 )  and  P ( f ! 5 ) •  Under  the  assumption  that  runs 
of  six  or  higher  need  not  be  evaluated  (because  their  contribution  to  the 
probability  of  failure  tends  to  zero),  no  more  diagrams  need  be  prepared.  Note 
for  checking  purposes  that  the  run  probabilities  of  i  jugs  being  outside  the 
crane  error  circle  sum  to  one.  The  terms  P(2),  etc.,  in  Equation  BIO  are 
found  as  usual  by  the  binomial  distribution, 


10  2  10-2 
P(2  out)  =  C  p  (1-p) 

2 

So  Equation  B10  becomes, 
P(failure)  cs  ( .2570) ( .05)  + 
—  0.0399  for  ten 


10  2  8 
=  C  (0.135)  (0.865)  =  0.2570. 

2 


( . 1070) ( .15)  +  ( . 02922 )( .2917)  +  (.00547) (.4583) 

jugs,  against  the  correct  value  of  0.0404 

in  Table  1.  (B14) 
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The  reduction  in  terms  made  in  Equation  BIO  is  argued  as  follows.  The 
probability  of  the  crane  mislaying  six  jugs  tends  to  zero,  e.g.,  P(6  out)  = 
0.000,  so  the  terms  pertaining  to  larger  jug  numbers  go  to  zero  and  the  sum 
of  those  higher  terms  does  not  contribute  to  P(failure).  In  fact,  the 
argument  is  false  as  computer  calculations  show.  The  practical  effect  of 
accepting  the  argument  is  to  mislead  one  into  thinking  that  long  juglines  are 
no  more  apt  to  fail  than  a  6- jug  line.  Actually,  probability  of  failure  is 
almost  linear  with  jug  numbers  up  to  30,  as  shown  in  Figure  6.  That  figure 
also  illustrates  that  the  rate  of  increase  of  probability  of  failure  is  greater 
when  p  =  0.2  than  when  p  =  0.1.  When  p  =  0.1,  the  30-jug  line  is  not  much  more 
apt  to  fail  than  a  6-jug  line  ( 7 %  vs.  1 %) ;  however,  when  p  =  0.2,  the  proba¬ 
bility  of  failure  goes  from  5$  for  six  jugs  to  25$  for  30  jugs.  The  upper 
slope  (p  =  0.2)  is  3.3  times  steeper  than  the  lower  slope  (p  =  0.1 ). 
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APPENDIX  C 
RUNS  OF  ARRANGEMENTS 
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In  Appendix  B  the  probability  of  failure  was  found  by  considering  every 
possible  arrangement  of  jugs.  Here  the  probability  of  failure  is  found  by 
considering  all  the  possible  arrangements  of  runs  (one  jug  or  consecutive 
jugs  out  of  place),  and  the  number  of  ways  each  of  these  arrangements  can 
occur.  A  probability  is  calculated  for  each  possible  arrangement  and 
summed  over  all  arrangements  to  give  the  overall  probability  of  failure. 


The  equations  of  this  method  are  further  explained.  The  crane  lays  each 
jug  in  or  out  of  the  allowable  error  circle  around  the  ideal  mark.  Like 
items  that  individually  have  a  two-state  existence  (in  or  out,  up  or  down, 
heads  or  tails)  and  are  placed  or  counted  together  have  a  probability  of  a 
number  of  each  state  arising.  This  is  named  a  Bernoulli  or  binomial  process. 
By  that  process  the  probability  P  of  laying  nQut  jugs  out  and  the  probability 
Q  of  laying  the  remaining  (n  -  nQut)  jugs  in  is  the  product  of  each 
independent  laying  event, 

n  ,  n— n  4. 

P(arrangement)  =  P  °U  Q  ,°U  where  Q  =  (1  -  P).  (Cl) 


This  equation  only  recognizes  the  count  of  so  many  jugs  out  and  the  comple¬ 
mentary  number  in.  It  does  not  recognize  the  order  of  the  outs  and  ins 
(permutations  or  arrangements  of  jugs). 


Several  different  arrangements  of  n  jugs  having  the  same  number  n  out 
are  possible.  Call  the  number  of  such  arrangements  n^.  All  of  these  arrange¬ 
ments  were  pictured  in  Appendix  A  for  small  jug  lines!  i.e.,  n  =  3*  4*  5. 
Drawing  more  was  cumbersome.  The  following  formula  was  found  which  gives  the 
number  of  arrangements  possible, 


nL  = 


rn  -  n  +  r 
out 


(C2) 


where  R  =  number  of  runs  in  the  arrangement  of  jugs.  Because  each  arrangement 
has  the  same  probability  of  arising,  we  multiply  Equations  Cl  and  C2  to  get 
the  probability  of  all  arrangements  of  R  runs  with  n  jugs  mislaid. 


P(arrangement) 


P) 


n-n 


out 


(C3) 


Only  some  arrangements  of  jugs  will  cause  partial  failure  to  detonate 
the  cloud  line.  The  bad  arrangements  (Figure  4)  are  ones  with  adjacent 
jugs  on  the  left  and  right  side  of  the  crane  error  circle.  Testing  and 
comparison  with  the  diagrams  of  Appendix  A  show  that  the  portion  of 
arrangements,  i.e.,  probability  of  failure,  which  will  give  partial  failure  is 

1  -  C(L  +  1 )/2**L]  (C4) 

for  one  run  of  length  L  or 

1  -  [1  -  Pa,)  ]  ( 1  -  P(L2)][1  -  P(L3)J  ...  (C5) 

for  multiple  runs  where  P(L1),  P(L2)  are  given  by  Equation  C3. 

The  reasoning  behind  Equation  C4  was  given  beside  Equation  B11. 
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The  probability  of  partial  failure  is  the  probability  of  the  dependent 
joint  events  of:  (l)  an  arrangement  arising  and  (2)  the  probability  of  failure 
with  that  arrangement.  This  statement  is  written  as 

P(failure  with  an  arrangement)  =  P(f,A)  =  P(fjA)  P(A),  (C6) 

where  A  is  any  arrangement  of  jugs  in  or  out  of  the  error  circle. 

The  total  probability  of  failure  is  the  sum  of  the  independent  arrangements' 
probabilities  of  failure, 

P(failure)  = 

For  comparison  with  the  method  of  Appendix  B,  Equation  B7>  the  3- jug 
line  problem  is  worked  again  in  Table  Cl.  For  comparison  with  the  method 
of  Appendix  A,  Equations  A9-A11,  a  5-jug  line  is  worked  again  in  Table  C2. 
These  calculations  use  the  probability  of  mislaying  one  jug  as  p  =  0.135. 

The  method  of  this  appendix  was  coded  in  FORTRAN  77  and  the  program, 
named  JUGRUN,  is  listed  in  Table  C3.  The  results  tables  were  generated 
by  the  program.  The  program  is  efficient  up  to  n  =  25  jugs  but  gets  very 
time  consuming  to  run  after  that.  Some  shorter,  faster  programs  based  on 
set  theory  (not  run  theory)  will  be  listed  in  Appendix  D. 

TABLE  Cl.  Probability  of  Failure,  Runs  of  Arrangements,  3-Jug  Line. 


P(failure  with  an  arrangement).  (C7) 


arrangements  of 

one  run  arrangements 

of  " 

1  or  2 

of  length  one 

runs  of 

length  2 

or  1 

R  =  1  +  + 

0  R  =  1 

+ 

0 

0 

L  =  1  +  0 

+  L  =  2 

0 

0 

+ 

0  + 

+ 

R  =  2 

0 

+ 

0 

L  =  1 

R  L 

no.  of  length 

P(f |L) 

nL 

P(L) 

runs  of  runs 

(Note  1) 

(Note  2) 

arrangement  of  one  run 
of  length  three 
R  *  1  0  0  0 

L  =  3 


P( failure ) 


(L+1) 


1 

1 

1 

2 


n  =  L 
out 

1 

2 

3 

1,1 


1  - 


n-n  +  1' 
out 


nLP(f iL)P(L) 


0 

3  0.1351 

0.8652= 

0.1010 

0 

1/4 

2  0.1352 

0.8651 = 

0.0158 

0.0079 

1/2 

1  0.1353 

0.865°= 

0.0025 

0.0012 

1-0)(l-0) 

1  0.1352 

0.865^= 

0.0158 

0 

run 

failure  for 

the  3-jug  line  is  0.0091  as 

0.0091 

shown  in  Table  1  and  earlier  appendices. 


Note  1 . 
Note  2. 


Probability  of  failure  given  this  run  length. 
Probability  of  getting  this  many  out  and  in. 


P(L)  = 


p  °ut  (i-p) 


I— n 


out 


Binomial  formula: 
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TABLE  C2.  Probability  of  Failure,  Runs  of  Arrangements,  5-Jug  Line. 


R  L 


no.  of 

runs 

length  P(f'L) 

of  runs  (Note  1) 

nL 

P(L)  P(failure) 

(Note  2) 

n  = 

out 

(L+1)  /n- 

1  ( 

'nout+ 

R 

) 

nLP(f|L)P(L) 

1 

1 

0 

5 

0.1351 

0.8654=  0.075578 

0 

1 

2 

1/4 

4 

0.1352 

0.8653=  0.011795 

0.011795 

1 

3 

1/2 

3 

0.1353 

0.8652=  0.001841 

0.002762 

1 

4 

11/16 

2 

0.1354 

0.865^  =  0.000287 

0.000395 

1 

5 

26/32 

1 

0.1355 

0.865°=  0.000045 

0.000037 

2 

1,1 

1 — ( 1 —0 ) ( 1 —0 ) 

6 

0.1352 

0.8653=  0.011795 

0 

2 

1,2 

1-(1-0)(l-l/4) 

3 

0.1353 

0.8652=  0.001841 

0.001381 

2 

1,3 

1-(1-1/4)(1-l/2) 

1 

0.1354 

0.8651®  0.000287 

0.000144 

2 

2,1 

1-(l-l/4) (1-1/2) 

3 

0.1353 

0.8652=  0.001841 

0.001381 

2 

2,2 

1— (1— 1/ 4) (1—1  /4) 

1 

0.1354 

0.8651-  0.000287 

0.000126 

2 

3,1 

1— (1—1/2) (1— 0) 

1 

0.1354 

0.8651*  0.000287 

0.000144 

3 

1,1,1 

1-(l-0)(l-0)(l-0) 

1 

0.1353 

0.8652=  0.001841 

0 

0.018165 

shown  in  Table  1:  0.0182 


Note  1.  Probability  of  failure  given  this  run  length. 

Note  2.  Probability  of  getting  this  many  out  and  in.  Binomial  formula: 

n  „  n-n 

P(L)  -  P  °Ut  (1-P)  °Ut 
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TABLE  C 3  .  JUGRUN  Program. 


PPO'.PAfi  T  ..  A L  C V  L  A 1  E  1  Ht.  P ft OB  A E?  ’  l  i  •  .  OP  FA  I  LURE  TO 
D£T  uNAT  E  DUE  T U  CFA’-E  F.SROF.  I  il !  -Pi..  a  .ED  JUGS  > 
i  n  P  L.  I  C I T  D  0  U  R  L  E  ft  u  I  C  I  0  N  :  A  -  H  G  -  Z  ' 

D  I  MENS  10  N  C  P  R  0  B  «  ? :  j  '  PR  0  p  t  ;•  0  1 
UR i t E ( 6  ,  :  1 

i  FORMAT < /'CALCULATION  OF  PROBABILITY  OF  FAIL.UPE  TC  DETONATE ' / 

* ' D •  j E  TO  CRaNE  ERRC  P  ‘Ml S P L , A C E 0  J 0 0 S )  / / 

’  WHEN  PROMPTED .  ENTER  N  THE  NUMBER  OF  JUGS  IN  A  ROW,.'/ 

=♦  'AND  P,  THE  P  P  0  8  A  P I  LIT  OF  AN  INDIVIDUAL  JUG  BEING  LAID'/ 

*  ’OUTSIDE  IT?  ERROR  CIRCLE  '/• 

* ' N  MUST  BE  A  POSITIVE  INTEGER  NOT  GREATER  THAN  30,'/ 

*  ’  AND  F  MUST  BE  GREATER  THAN  f*  AND  LESS  THAN  1  '/' 

* 'ENTER  ZEROES  FOR  N  AND  P  WHEN  NO  MORE  DATA  ARE  TO  BE  ENTERED  '/) 

£  READ  I  S ,  •*  )  N  ,  P 

IF  IN  EG  0 ) GO  TO  £! 

COMPUTE  CONDITIONAL  PROBABILITY  OF  FAILURE  GIVEN 
I  CONSECUTIVE  JUGS  ARE  OUTSIDE  THEIR  ERROR  CIRCLES, 

AND  PROBABILITY  OF  I  JUGS  BEING  PLACED  OUTSIDE 
THEIR  ERROR  CIRCLES,  FOR  1=1.0,  .  N 

0=1  o-p 
DO  3  1=1 , N 

CPRQBi.  t  )  -  1  0- (  1  +  1  .  0  ) /£**  I 

3  PROP  i  I  1  v-p  * *  I+Q  +  MN-!  i 
PRFAIL=0  0 

CALCULATIONS  FOR  1  RUN 
DO  4  NOUT = ! . N 
NFREG  =  NC0M9( N  +  1 -NOUT ,  1  ) 

LRUN1 =NOUT 

PROBC  =  CF'ROB  ( LRUN1  ) 

4  RRFAIL=PRFAIL+NFRE9+PRU8C*RR0&: NOUT ) 

CALCULATIONS  FOR  £  RUNS  (MUST  HAVE  AT  LEAST  3  JUGS) 

IF- N  LT  3) GO  TO  19 

DO  S  N 0 U T  =  £ , N - 1 

NFREQ  =  NC0f1B  (  N+  1  -NOUT  ,  2  ) 

DO  S  L  RUN  1=1. NOUT - 1 
LRUN£=NOUT- LRUN1 

P R 0 B C  = 1  . 0-i  1  . O-CPROBCLRUNI  >  )  +  M  O-CPROB ( LRUN£ >  i 

5  P  R F  A  I L  =  P R F  A  l  L  +  M F R E 0 F  R 0 B C  t-PROE  f  N0< 'T  ) 

calculations  for  3  runs  t must  have  at  least  5  jugs > 

IF' N  LT . £ ) GO  TO  1 9 
DO  6  N 1 0 T  =  3  N— >* 

NFREO  =  NCOMB ( N+ i -NOUT  .  3 ) 

DO  »>  s.SUN1  -1  ,  NOUT -2 

DC  (■>  LRUN£=l  N0UT-1-LRUN1 

L  P  U  N  3  =  N  0  L*  T  -  L  R  U  N  1  -  L  R  U  N  c 

PRG8C  *  1  0-'.  1  0  •  •  C  P  R  0  B  (  w  R  U  N 1  )  )  *(  1  0 -CPROB  (  LRUN2  )  >  M  1  .  O-CPROE  (  LRUN3  )  1 
C  P R  F a i L  =  F  R  F  A  I l N P  R E 0  +  F  ft 0 b C  +  P ROD ( N 0 1  j T ) 
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TABLE  C3.  (continued) 


CALCULATIONS  FOR  4  RUNS  i  MU  3  i  HAVE  AT  LEAST  7  JUGS  • 

IF i N  LT  GO  TO  13 
DO  7  NOUT-4 , N-3 
NFREQ=NC0M8( N+1  -NO'JT  ,  4  ) 

DO  7  LRUN1 =1 , NOUT-3 

DO  7  LRUN2-1  , N0UT-2-LP UN  1 

DO  7  LRUN3=t  ,  NO'JT  - 1 -LRUN  1 -LRUN2 

L  P  U N 4  =  NO U T — L R UN  t -LRUN2-LRUN3 

r  R  0  &  C  —  i  0  -  (  1  ■  0  -  C  r  R  0  B  (  L  S  U  N 1  5  )  *  '.  1  .  O-CPROB  (  LRUNc'  '  )  * 1  1  O-CPRGBl  LRUN3  )  ) 
+  *  <  1  O-CPROB i LRUN4  )  ) 

7  FRF  A  IL  =  PRFAIL  +  NFREQ.*PR0BC*PR08(N0UT  ) 

CALCULATIONS  FOR  5  PUNS  (HOST  HAVE  AT  LEAST  9  JUGS  5 

t  F • H  LT.9)00  TO  19 

L1 0  3  N  OUT  —  5  ,  N  - 

NFPtO-NCOMB(N+1-NOUT ,55 

DO  8  LRUN1 - 1 , NOUT-4 

DO  8  LRUN2=1 , N0UT-3-LRUN1 

DO  8  LRUN3=1 , N0UT-2-LRUN1 -LRUN2 

DO  8  LRUN4=1 , NOUT- 1 -LRUN1 -LRUN2-LRUN3 

LRUN5=N0UT-LRUN1 -LRUN2-LRUN3-LRUN4 

PROBC- 1  . 0-(  1  0- CPROB(LRUNI)  )  +  ( 1  . O-CPROB ( LRUN2 5  )  +  ( 1  . O-CPROB ( LRUN3 ) 5 
**(  I  .  O-CPROB ( LRUN4 ) >*( 1  . O-CPROB ( LRUN5 ) ) 

8  PRFAIL=PRFAIL+NFREQ*PROBC*PROB(NOUT ) 

CALCULATIONS  FOR  6  RUNS  (MUST  HAVE  AT  LEAST  11  JUGS  5 

I  F ( N  LT .  1  1  )G0  TO  19 

DO  9  NOUT  =  6 , N-5 

NFREQ=NC0MB(N+1 -NOUT , 6 5 

DO  9  L.RUN1  =1  ,  NOUT-5 

DO  9  LRUN2=1 , N0UT-4-LRUN 1 

DO  9  LRUN3=1 , N0UT-3-LRUN1 -LRUN2 

DO  9  LRUN4=1 , N0UT-2-LRUN1 -LRUN2-LRUN3 

DO  9  LRUN5= 1 , NOyT- 1 -LRUN 1 -LRUN2-LRUN3-LRUN4 

LRUN6=N0UT-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5 

P  R  0  B  C  =  1  0 -  l  1  O-CPROB ( LRUN 1  ) )*<  1  O-CPROB { LRUN2 )  )*( 1  O-CPROB ( LRUN3 ) ) 
**(  1  .  O-CPROB  <  LR.UN4  )  )*(  1  .  O-CPROB  (  LRUN5  )  )*(  1  O-CPROB  (  LR  UNA  )  5 

9  PRFAIL=PRFAIL*NFREQ*PROBC*PROB(NOUT 5 

CALCULAT I ONS  FOR  7  RUNS  (MUST  HAVE  AT  LEAST  13  JUGS) 

IF . N . LT . 1 3) GO  TO  19 
DO  10  N0UT=7 , N-G 
NFREQ=NC0MB(N+1 -NOUT . 7) 

DO  1 0  LRUN 1 -1 , NCUT-S 
DO  10  LRUN2- 1 , N0UT-5-LRUN ■ 

DO  10  LRUN3=1  ,  N0UT-4-L.RUN  1 -LRUN2 

DO  10  LRUN4=1 , N0UT-3-LRUN1 -LRUN2-LRUN3 

DO  10  LRUN5= 1 , NOUT- 2- LRUN 1 -LRUN2-LRUN3-LRUN4 

DO  10  LRUN6=1  ,N0UT-1-LRUN1-LRUN2-LRUN3-LRUN4-LRUf:t 

LRUN7=N0UT-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6 

PROBC=1  0-  (  1  O-CPROB  (  LRUN  1  )  )  *(  1  O-CPROB  (  L  PUN3  5  i  •*  <  ■  "-CFROP  i  L.RUN3  )  > 
t*(  1  O-CPROB  <  LRUN4  5  ;  +  (  1  .  O-CPROBi  LR'JNS)  )*<  1  O-CFROE  •  t.PUN6  )  >  *• 

Ml  0-CPROB(LRUN7 >  ) 

1  0  PRFAIL=PRFAIL  +  NFREQt‘P:>OBC  +  pROB(NOUT) 


43 


TABLE  C3.  (continued) 

C  CALCULATIONS  FOR  8  RUNS  <  MUST  HA'-'E  AT  LEAST  IS  JUGS' 

IFlN.LT. IS) GO  TO  19 
DO  11  NOUT  =  8 N-7 
NFRE0  =  NC0MB(N+1-N0UT,8.» 

DO  1 1  LRUN1 -1 , NOUT-7 
DO  11  LRUN£=1 , NOUT-6-LRUN 1 
DO  11  LRUN3=1  , NQUT-5-L RUN  1 -LRUN2 
DO  11  LRUN4  =1,N OUT-4- I..R UNI -LRUNS-LRUN3 
DO  11  LRUNE-1  .  NOUT-2 -LR.UM1  -  LFUMC-  L  RUNT -LPUN4 
DO  11  LRUN6=1 . N0UT-2-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5 
DO  1  1  LRUN7= 1  , NOUT- 1 -LRUN1 -LRUN2-LRUN3-LRUN4-LRUNS-LPUN6 
LRUN8  =  N0UT-LRUN1  -LRUN2-LRUN3-L.RUN4-LRUN5-LRUNE-LRUN7 
PROBC* 1  0 - ( 1 . 0-CPR0B1LRUN1 > )*( !  0-CPRGB ( LRUN2 ' )*M  . O-OPROB t LRUN3 ) > 
*+(  l  O-CPROB  (  LRUN4  )  i  1  O-CPROB ■ LPUN5 >  '*(1  O-CPROB <  LRUN6 >  >* 

*(1  O-CPROB i LRUN7  )  )  *( 1  O-CPROB ( LRUN8  )  ) 

1 1  PRFAIL=PRFAIL+NFREQ*PROBC*PROS(NOUT> 

C  CALCULATIONS  FOP  9  RUNS  (MUST  HAVE  AT  LEAST  17  JUGS' 

IF1N.LT. 1 7  5  GO  TO  19 
DO  12  NOUT=9 , N-8 
NFREQ=NCOMB( N+1 -NOUT , 9 ) 

DO  1  £  LRUN1 =1 , NOUT-3 

DO  1  £  LRUN2=1  , N0UT-7-LRUN 1 

DO  12  LRUN3=1 , NOUT-6-LRUN1 -LRUN2 

DO  12  LRUN4=1 , N0UT-5-LRUN 1 -LRUN2-LRUN3 

DO  12  LRUN5=1 , N0UT-4-LRUN 1 -LRUN2-LRUN3-LRUN4 

DO  12  LRUN6=1 , N0UT-3-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN3 

DO  12  LRUN7=1 , N0UT-2-LRUN1 -LRUN2-LRUN3-LRUN4-LRUNS-LRUN6 

DO  12  L  R  U  N  8  =  1  , NOUT- 1  - LRUN 1 -LRUN2-L RUN3-LRUN4- L RUNS-LRUN6- LRUNT 

LRUN9=N0UT-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUNfc~LR'JN7--LRUN« 

PROBC  =  1  0-  (  1  .  O-CPROB  <  LR'JNl  ))♦< 1  .  O-CPROB  (  LRUN2  )  )  *  <  1  O-CPROB  (  LRUN3  )  ' 
**(1  O-CPROB i LRUN4 ) )*< 1  . 0 -CPROE ( LRUN5 ) >*l  1  . O-CPROBI LRUN6 ) ) + 

* <  1  . 0-CPROB( LRUN7  )  ) * ( 1  O-CPROB < LRUNB ))*( 1  O-CPROB < LRUN9  )  ) 

1  £  PRFAIL-PRFAIL+NFREQ*PRGBC*PRGB( NOUT  ) 

C  CALCULATIONS  FOR  10  RUNS  (MUST  HAVE  AT  LEAST  19  JUGS) 

I F ( N  LT .  1 9  >GO  TO  1 9 
DO  1 3  NOUT- 1 0 , N-9 
NFREQ  =  NCOMB(  N+1  -NOUT..  1  0  ) 

DO  1 3  LRUN 1 =1 , NOUT-9 

DO  13  LRUNE=1 , N0UT-8-LRUN 1 

DO  13  LRUN3= I , N0UT-7-LRUN1 -LRUN2 

DO  13  LRUN4=1  ,  NOUT-6-LRUN1 -LRUN2--LRUN3 

DO  13  LRUN5=I , N0UT-5-LRUN1 -LRUN2-LRUN3-LRUN4 

DO  13  LRUNC=1 , N0UT-4-LRUN1 -LPUN2-LRUN3-LPUN4-LRUNS 

DO  13  LRUN7=1 , NOUT -3-LRUN 1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN& 

DO  13  LRUNS= 1 , NOU 7 -£- LRUN 1-LPUN2-LRUN3-LRUN4- LRUN 3-LRUN6 -LRUN? 

DO  13  LRUN9  -  1  ..  NOUT- 1 -LRUN  1  -  LRUN2-LRUN3-LRUN4-LRUNE-LRUN& -LRUN7- 
+LRUN8 

LRUN 1 0=N0UT-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7-LRUN8-LRUN9 
PROBC=  1  O-i  1  O-CFROBi  LRUN1  '  »*(  1  .  O-CPROB  (  LRUN2  )  )*(  1  O-CPROB  (  LRUN3  )  .' 
**( 1  O-CPROB I LRUN4 ) )  *  (  1  . O-CFROB ( LRUN5 > )  +  (  1  O-CPROB  <  LRUN6 ) )* 

O-CPROB  (  LRUN7  ))*(  1  O-CPROB  (  LRUN8  ))*<  1  O-CPROB ( LRUN9 ) M 
*  (  1  0  -CPROB  •  LPIJN1  0  >  ) 

1  3  PRF AIL  =  PRFA l  L+NFREO  +  PROBC +  PROE ( NOUT ) 
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TABLE  C3.  (continued) 


L  l-ALLULA  i  i UNS  FOR  11  RUNS  (MUST  HAVE  AT  LEAST  £1  JUGS) 

I F ( N . LT . £  1  ) GO  TO  19 
DO  1  4  NOUT  = 1  1  . N-1 0 
NFREQ=NCOMB(N+1 -NOUT, 1 1  ) 

DO  1  A  LRUN1 =1  , NOUT- 1 0 

DO  14  LRUN2=1 , NOUT-9-LRUN 1 

DO  14  LRUN3=1 , N0UT-8-LRUN1 -LRUN£ 

DO  14  LRUN4=1 , N0UT-7-LRUN1 -LRUN2-LRUN3 
DO  14  LRUNS=1 , NOUT-6-LRUN1 -LRUN2-LRUN3-LRUN4 
DO  14  LRUN6=1 , N0UT-5-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5 
DO  14  LRUN7=1 , N0UT-4-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6 
DO  14  LRUN8=1 ,N0UT-3-LRUN1-LRUN£-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7 
DO  14  LRUN9= 1 , NOUT -2-LRUN 1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7- 
*LRUN8 

DO  14  LRUN1 0=1 ,N0UT-1-LRUN1-LRUN2-LRUN3-LRUN4-LRUN5-LRUN6~LRUN7- 
+LRUN8-LRUN9 

LRUN1  1  =NOUT-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7-LRUN8-LRUN9- 

♦  LRUN 1 0 

PROBC= 1  . 0-( 1  . O-CPROB ( LRUN1 ) )*(  1  . O-CPROB ( LRUN£  )  )*(  1  . O-CPROB ( LRUN3  )  ) 
**( 1  . O-CPROB (LRUN4)  )*<  1  . O-CPROB ( LRUN5 )  )*( 1  . O-CPROB f LRUN6  )  )* 

*(  1  . O-CPROB ( LRUN7  )  )  *( 1  . O-CPROB ( LRUN8  )  >*( 1  . O-CPROB ( LRUN9  )  )* 

*( 1  . O-CPROB (LRUN1 0  )  )*( 1  . O-CPROB ( LRUN 1 1  )  ) 

14  PRFAIL=PRFAIL+NFREQ*PROBC*PROB(NOUT) 

C  CALCULATIONS  FOR  1£  RUNS  (MUST  HAVE  AT  LEAST  £3  JUGS) 

I F ( N . LT . £3 ) GO  TO  19 
DO  15  NOUT  =  1  £ , N-  1  1 
NFREQ=NCOMB ( N+ 1 -NOUT ,12) 

DO  15  LRUN 1 = 1 , NOUT- 1 1 
DO  15  LRUN2=1 , NOUT- 1 0-LRUN 1 
DO  15  LRUN3=1 , NOUT-9-LRUN! -LRUN2 
DO  15  LRUN4=1 , N0UT-8-LRUN 1 -LRUN2-LRUN3 
DO  15  LRUN5=1 , N0UT-7-LRUN1 -LRUN2-LRUN3-LRUN4 
DO  15  LRUN6=1 , NOUT-6-LRUN 1 -LRUN2-LRUN3-LRUN4-LRUN5 
DO  15  LRUN 7=1 , N0UT-5-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6 
DO  15  LRUN8=1 ,N0UT-4-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7 
DO  15  LRUN9=1 ,N0UT-3-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7- 
*LRUN8 

DO  15  LRUN 1 0=1 ,NOUT-£-LRUN1 -LRUN£-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7- 
♦LRUN8-LRUN9 

DO  15  LRUN1 1=1 , NOUT-1-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7- 

♦  LRUN8-LRUN9-I.RUN1  0 

LRUN 1 £  =  NOUr -LRUN I -LRUN2 -LRUN3 -LRUN4-LRUN5 -LRUN6 -LRUN7 -LRUN8  - -LPUN9 - 
♦LRUN 1 0-LRUN 1 1 

PROBC  =  1  . 0- ( 1  . O-CPROB ( LRUN 1  ) )♦( 1  . O-CPROB ( LRUN2 ) )*(  1  . O-CPROB ( LRUN3 ) ) 
**( 1  . O-CPROB ( LRUN4  )  )*( 1  . O-CPROB ( LRUN5  )  )*( 1  O-CPROB ( LRUN6  )  )* 

*< 1  O-CPROB ( LRUN? ) )*( 1  . 0  - CPROB ( LRUN8 ) )*( 1  . O-CPROB ( LRUN9 ) >* 

♦  ( 1  . O-CPROB ( LRUN 1 0 ) )*( 1  . 0 -CPROB ( LRUN  1  1  ) )*(  1  O-CPROB ( LRUN 1  £ ) ) 

15  PRFAIL=PRFAIL+NFREQ*PROBC*PRCB(NQUT) 
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TABLE  C3.  (continued) 


0  CALCULATION-*  FOR  13  RUN'S  MUST  KAv'E  AT  LEAST  2S  JUGS1 
IFfN.LT  25)00  TO  IT 
DC.  16  N 0 U T  =  1  3  .  N  -  1  it 
NFRE(5  =  NCOMP  (  N+?  -NOUT  1  3  '• 

DO  IS  LRUNI = I . NOUT- 1 2 

DO  16  LRUNE=1 , NOUT-1 ! -LRUN1 

DO  16  LRUN3=  1  ;  NOUT- 1  0- LP'JN  I -LPUNZ 

DO  16  I.RUN4  -  1  NOUT-9-LRUN1 -LPUN2-  LRUN3 

DO  16  LPUNE  -  1  .  NO'JT-8-LRUNI -LPUN£-;.R«JN.3-LPUN4 

DO  16  L RUN 6-1  NOVT -7- LR'JN  '  -LRUN2-LRUN2-LRUN4 -!  RUN" 

DO  16  '  RUN? "  1  -NOU  r  6 - L R U N 1  - L  F  UN  £ - LR 1  j  N  3 - L  R  U  N  4 - i  RUNT-LRUN6 
DO  16  LRUNS-1  ,  N0UT-5-LRUN1 -LRUNP-L  RUN3-LPUN4-LRUNS-LRUN&-'.RUN7 
DO  16  L  R  U  N  ?  »  1  .  N  G  i  j  >  -  4  - 1  RUN  i  -LRf.'Nt-LPUN  i-LRUN4-i  P>JN~-LRUN6  -l.  RUNT- 

*  L.  K  U  N  3 

DO  16  LRUNI 0=1 . N0UT-3- LRUN1 -LRUN2-LR UN 3-LRUN 4~ LRUN5- LB UN G -LRUN7- 
♦LRUN8-LRUN9 

DO  16  LRUN1  1=1  , N0UT-£~LRUN1-LRUN2-LRUN3-LRUN4-LPUNS-LPUNb-LRUN7- 

*  LRUNS-1.  RUN9-LPUN1  0 

DO  16  LRUN12-1  ,  NOU  T-1 -LRUN1 -LRUN8-LRUN3-LRUN4-LRUNE.-LRUN6-LRUN7- 
♦LRUN8-LRUN9-LRUN1 0-LRUN1  1 

LRUN13=N0UT-LRUN1-LRUN£-LRUN3~LRUN4-LRUN£-LRUNG-LRUN7-LRUN8-LRUN9- 
+LRUN1 0-LRUNl 1 -LRUN1 2 

PRCBC=  t  .  0-f  1  . O-CPROB- LRUN1  ) )*(  1  . O-CPROB ( LPUN2  )  )*( 1  . O-CPROB ( LRUN3 ) J 
*+( 1  . O-CPROB l LPUN4  )  )*(  1  . O-CPROB  <  LRUN5  )  )*(  1  . O-CPROB ( LRUN6  )  )* 

*( 1  . O-CPROB ( LRUN7  )  )  M  1  . O-CPROB ( LRUN8  )  )*< 1  O-CPROB ( uRUN9  )  )  * 

•M 1  O-CPROB ( LRUN 1 0)  )*( J  O-CPROB ( LRUNI 1  )  )*( 1  0-CPR05 ( LfiUN 1  £  )  ) * 

Ml.  O-CPROB  (  LRUN1  3)  ) 

1  6  PRFAIL  =  PRFAIL  +  NFREQ*PROBC*PROP<NO»'T) 

C  CALCULATIONS  FOR  14  RUNS  (MUST  HAVE  AT  LEAST  £7  JUGS) 

IF(N  LT  £7) GO  TO  19 
DO  1 7  NOUT  =  1  4  ,  N-  1  3 
NFREQ  =  NCOMB< N+1 -NOUT,  14) 

DO  1 7  LRUN1 =1 , NOUT-1 3 

DO  17  LRUN2=1 ,N0UT-12~LRUN1 

DO  17  LRUN3=1  ,  NOUT-1  1 -I.RUN1 -LRUN2 

DO  17  LRUN4=1 , NOUT-1 O-LRUN 1 -LRUN2-LRUN3 

DO  17  LRUN5=1  , NOUT- 9- L RUN  1 -LRUN£-LRUN3-LRUN4 

DO  17  LRUN6=1 , N0UT-8-LRUN 1 -LRUN2-LRUN3-LRUN4-LRUN5 

DO  17  LRUN 7=1 , NOUT -7 -LRUN 1 -LRUN2-LRUN3-LRUN4-LRUNE-LRUN6 

DO  17  LR'JN 8  =  1  ,  NOUT- 6 -LRUN  1  -LPUN2-l.RUN3-LRUN4-LRUN5-LRUN6-L.RUN7 

DO  1  7  L R U N 9  = 1  , N0UT-5-LRUN 1 -LRUN2-LRUN3-LRUN4-LRUNS-LRUN6-LRUN7- 

*  L.  R  i;  N  8 

DO  17  LRUN I  0  =  1  . N0UT-4-LRUN1 -LPUN2-LRUN3-LPUN4-LRUNB -LRUNt-LRUN? - 
*LRUN3 -LRUN9 

DO  1  7  LRUN  1  1=1,  N0UT-3-LRUN 1  -LPUN2  -LRUN3-LRUN4 - I.RUN5-LRUN6  - LRUN7- 

*  L  R  U N 8 - 1.  PUNS  -LRUN  1  0 

DO  17  LRUM3-1  ,NOUT-£-LPUN1-LRUN£LRUN3'LRUN4-LRUNS-LPUN6-LRUN7- 

*  LRUN8-LPUN9-LRUN1 O-lRUN  I  1 

DO  17  [.RUN  1  3=  1  .  NOUT-1  -LRUN1 -LRUNE-LRUN3-LRUN4-LRUN5-LRUN6--LRUN7- 
*LRUN8-LRUN9-LPWNt 0 -LRUN 1 1 -LRUN 1 2 

LRUNi  ^ - N G U T  -  LRUN  1  - LRUN2-LRUN3 -LRUN4-LRUNE- LRUN6-LRUN7-I.RUN8-LRUN9- 
v-LRUNl  O-LRUN  ■  1  -LRUNI 2-LRUN 1 3 

PROBE  =1  0-M  0- CPFOP  '•  LRUNI  )  )  M  I  0-CPRGB  f  LRUN£  )  )  M  1  .  0 -CPRCB  <  LRUN3  )  ) 
f  h  1  O-'.PROB  (  LPUN4  )  i  *  l  1  .  O-CPROB  <  LRUNS  /  )  M  1  O-CPROB  :  LRUN 6  1  )  * 

*  '  !  O-CPROB  ■'  LRUN  7  ))*d.  0  -CPROB  f  LRUN8  )  >  +  ;  1  O-CPROB  (  LRUN9  )  >  * 

<  1  o  -i  PR  OBi  l_  R  U  N  «  C  »  >*n  O-CPROF:  f  lRI.N1  1  >  )*(  1  .  O-CPROB  (  LRUN  1  2  »  )  * 
t  i*  0  Cf'ROB  ■:  i.RUN  13i)t  (  i  O-CPROB  (  lRUN  14?) 

’  7  P  P  F  a  I -PST  I L  +  N  F  R  E  Q  *  P  P  0  B  C  *  PR  OB  •  N U  U  T  ) 
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TABLE  C3.  (continued) 


CAL  CU  i  4TION5  FOR  1  £  RUNF  (  M»KT  MAWF  AT  i  FACT  JMr.c.j 
i  !-  t  N  .  L  T  c 9  )  G  U  T  U  I  9 
DO  18  NOUT  =  1  5  ,  N-  M 
NFREQ=NC0MB(N+1 -NOUT, 15) 

DO  1 8  LRUN1 =1 , NOUT-1 4 
DO  18  LRUN2=1 , NOUT-1 3-LRUN1 
DO  18  LRUN3=1 , NOUT-1 2-LRUN1 -LRUN2 
DO  18  LRUN4= 1 , NOUT-1 1-LRUN1-LRUN£-LRUN3 
DO  18  LRUN5=t , NOUT-1 0-LRUN1 -LRUN2-LRUN3-LRUN4 
DO  18  LRUN6=1 , NOUT-9-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5 
DO  18  LRUN7=1 , N0UT-8-LRUN1 -LRUN2-LRUN3-LPUN4-LRUN5-LRUN6 
DO  18  LRUN8=  1  , N0UT-7-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7 
DO  18  LRUN9=1 , N0UT-6-LRUN1 -LRUN2-LRUN3-LRUN4-LPUN5-LRUN6-LRUN7- 
*LRUN8 

DO  18  LRUN 10=1 . N0UT-5-LRUN1 -LRUNS-LRUN3-LRUN4 -LRUN5-LRUN6-LRUN7- 
♦LRUN8-LRUN9 

DO  18  LRUN 1 1=1 , N0UT-4-LRUN1 -LRUN2-LRUN3-LRUN4- LRUN5-LRUN6-LRUN7- 
♦LRUN8-LRUN9-LPUN1 0 

DO  18  LRUN 1 £= 1 , NOUT-3-LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7- 
*LRUN8~LPUN9~LRUN1 0-LRUN1  1 

DO  18  LRUN 1 3=1 , N0UT-2-LRUN 1 -LRUN2-LRUN3-LRUN4-LRUN5-L RUN6-LRUN7- 
*LRUN8-L.RUN9-LRUN1  0-LRUN1  1  -LRUN  1  £ 

DO  18  LRUN 14=1 , NOUT-1 -LRUN1 -LRUN2-LRUN3-LRUN4-LRUN5-LRUN6-LRUN7- 
♦LRUN8-LRUN9-LRUN1 0-LRUN1 1 -LRUN 1 2-LRUN 1 3 
LRUN15=N0UT-LRUN1-LRUN2-LRUN3-LRUN4-LRUN5-LRUN6~LRUN7-LRUN8-LRUN9- 
+  L  RUN  1 0-L  RUN  1 1 -LRUN 1 2-LRUN  1  3-LRUN 1 4 
PRQBC=1  . 0-( 1  O-CPROB ( LRUN 1  ) >  *  M  . 0- CP R OB  (  LRUN 2 » )  *  (  1  O-CPROB ( LRUN3  )  ) 
**( 1  O-CPROB (LRUN4)  )*<  I  . O-CPROB ( LRUN5 >  )*<  1  O-CPROB' LRUN 6 )  )* 

♦  (1  O-CPROB (LRUN7) ) * ( 1  . O-CPROB  f LRUN8 ) >*< 1  O-CPROB < LRUN9 )  >* 

*< 1  . O-CPROB (LRUN1 0) )*( i  . O-CPROB ( LRUN 1 1  ) )*M  . O-C PROP ( L RUN  1  £ ) )* 

*( 1  Q-CPROB(LRUN1 3) )*< 1  O-CPROB (LRUN14) )*<  1  O-CPROB! LRUN1 5 ) ) 

1 8  PRFAIL=PRFAIL+NFREQ*PROBC*PROBiN0UT ) 

19  URITE(6,£0)N,P .PRFAIL 

20  FORMAT ( / ' N  =',T18,I2/ 

+  ' P  = '  , T 1 4 , F  6 . 4// 

* ' PROBABILITY ' / 

* ' OF  FA ILURF  =  '  , F6  4/  ) 

GO  TO  2 

£1  STOP 
END 

FUNCTION  NCOMB (  N  ,  M  ) 

Ml =N-M 

IF ( M . LE . Ml  1  THEN 
MM-M 
ELSE 
MM-M  1 
END  IF 
NCOMB= 1 

IF (MM . NE . 0 ) THEN 
DO  1  1=1, MM 

1  NCOMB=NCOMB*(N+1 -I  )  f\ 

END  IF 
END 
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APPENDIX  D 
THE  BAD  EVENT 
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The  JUGFAE  probability  of  failure  due  to  improper  crane  placement 
can  be  computed  without  using  conditional  probability  or  run  theory. 

Two  easily  programmable  methods  for  doing  exact  calculations  are  given 
as  well  as  an  approximate  method  which  would  be  useful  for  field  work 
when  computing  facilities  are  not  available. 

The  first  method,  discussed  in  Section  D.1,  runs  into  some  numerical 
stability  problems  for  very  large  n  and  large  p  (close  to  1 )  but  these 
cases  are  operationally  unimportant.  For  the  value  of  p  (probability 
of  misplacing  one  jug)  that  is  considered  in  this  work  (p  =  0.135) 
the  method  works  for  any  n.  The  program  listing  of  this  method,  called 
JUGCOM,  is  in  Table  D1 . 

The  second  method  avoids  the  stability  problem,  although  for  a  given 
p  all  failure  probabilities  for  1,  2,  ...»  n-1  must  be  computed  in  order 
to  arrive  at  the  n-th  probability.  This  recursion  is  only  a  minor 
problem  in  terms  of  computer  time.  By  comparison,  the  run  theory  JUGRUN 
program  is  only  valid  for  n  <  31  and  takes  much  more  time.  The  program 
listing  of  this  method,  called  JUGREC,  is  in  Table  D2. 

D.1  Exact  calculation,  JUGCOM.  We  have  a  failure  if  one  jug  is  placed 
too  far  to  the  left  and  the  next  one  (to  the  right  of  the  first  one)  is 
placed  too  far  to  the  right.  Thus,  we  have  a  failure  in  a  line  of  n  jugs 
if  the  following  occurs: 

(jug  1  too  far  left  and  jug  2  too  far  right)  or 
(jug  2  too  far  left  and  jug  3  too  far  right)  or 


(jug  n-2  too  far  left  and  jug  n-1  too  far  right)  or 
(jug  n-1  too  far  left  and  jug  n  too  far  right). 

The  placement  of  each  jug  is  assumed  to  be  independent  of  the  placement 
of  the  other  jugs. 


The  probability  of  failure  for  an  n-jug  line  may  be  written  as 
P(failure)  =  P(E1  U  E2  U  ...  U  En_2  U  En_1), 


(Dl) 


where  E^  -  the  event:  jug  i  too  far  left  and  jug  i+1  too  far  right, 
for  i  =  1,  2,  ...,  n-1  and  " U"  means  "or". 

Note  that  events  E^  and  E2  are  defined  quite  differently  in  Appendix  A. 

Equation  3  indicated  the  expansion  of  Equation  2  for  two  sets;  for  n  sets, 

n-I  n-2 

P(E1U  E2  U  ...  UEn_2  U  En-1 )  =y^P(Ei)  Z  (— 1 ) ^ ( j— way  intersections),  (p2) 

j ^1 
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P(E^U  %2  U  •  •  •  UEn-2  U  En_i )  =  individual  probabilities 

-  1-way  intersections 
+  2-way  intersections 


n  2 

+  (-1)  [(n-2)-way  intersections], 


(D3) 


=  [P(E1 )  +  P(E2)  +  ...  +  P(En_2)  +  P(En_.,)] 

-  [p(e1  n  E2)  -  p(e1  n  e3)  -  ...  -  p(En_2  n  En_1)] 
+  [P(E1  n  e2  HE3)  +  ...  +  P(En_3  n  En_2  nE^)] 


+  (-i)n-2[p(E1  n  e2  n  ...  n  e„_2  n  En_1 ) J , 

where  "U"  means  "or"  and  "O"  means  "and". 

Not  all  the  E^'s  are  independent  events.  In  particular,  E^  and  E^+^  are 
mutually  exclusive  events  because  E.  and  E.  .  cannot  both  happen.  (Jug  i+1 
cannot  be  both  too  far  to  the  right  and  too  far  to  the  left. )  However,  all 
other  subsets  of  events  are  independent,  e.g.,  E^  and  Ei+2  are  independent, 
E^  and  E^+^  are  independent,  etc.  In  general,  all  intersections  of  E^  and 
Ef+i  events  will  be  zero  and  all  other  intersections  will  be  nonzero.  All 
cases  can  be  generalized  in  an  expression  for  probability  of  failure  P(f), 

in/2] 


(E>4) 


P(f)  = 


■2>‘T;V- 


(D5) 


where  1=1 

n  =  number  of  jugs  in  a  row, 

4p/2j  =  the  integer  part  of  n/2, 

p  =  probability  of  an  individual  jug  being  outside  error  circle, 

p/2  =  probability  of  an  individual  jug  being  too  far  left 

=  probability  of  an  individual  jug  being  too  far  right, 

p'  =  (p/2)^  -  probability  of  one  jug  being  too  far  left  and  the  next 
jug  being  too  far  right  (by  independence  of  placements  of  jugs) 

=  probability  of  a  particular  pair  of  consecutive  jugs  causing  failure. 

For  computing  purposes,  rename  the  quantity  after  the  summation  sign  T. , 

Ln/2J 


P(0  =  £ 
i=l 


T. 

i 


(D6) 


A  ?  ^cursive  formula  can  be  obtained  by  substituting  i-1  for  i  in  T.  and 
finding  the  ratio  Ti/ii+1.  Letting  TQ  =  -1 ,  we  have  1 


-(n-2i+1 )(n-2i+2) 


T.  = 

l 


(P’)(Ti_l)»  i  =  2 .  [n/2J. 


i(n-i+1 ) 

Equations  D6  and  D7  form  the  basis  of  the  program  JUGCOM. 


(f  7) 
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As  an  example,  we  will  handle  a  5- jug  line  by  expanding  both  Equation  D4 
and  Equation  D5. 


a.  Five- jug  line  by  Equation  D4« 

P(f)  =  P(E  )  +  P(E  )  +  P(E  )  +  P(E  ) 

1  2  3.  4  -^0  -rO 

-  P(E JZrfT)  -  P(E  HE  )  -  P(E  HE)-  P^PftT)  -  P(E  HE  )  - 
—  1  2  '13  14  " 


2  3 


2  4 


3  4 


p(e  )  +  p(e  n^rtE^  +  p(Ep-n ^eT )  +  p(ej>s-i^ ) 

2  3  2  4  ' — 1  3  4  •'"2  3  4 


-  p(e  n  e 


0 

E  )•  (D8) 

3  4 

Since  adjacent  events  E^  and  Ej+i  are  mutually  exclusive,  the  intersections 

of  such  events  are  the  null  set  and  the  probability  of  the  null  set  is  zero. 

Then  Equation  D8  becomes 

P(f)  =  P(E  )  +  P(E  )  +  P(E  )  +  P(E  )  -  P(E  E  )  -  P(E  E  )  -  P(E  E  ),  (D9) 

1  2  3  4  13  14  24 

but  PCE^)  =  p'f  and  P(independent  events)  is  the  product  of  the  individual 


probabilities,  so  Equation  D9  becomes 

2 

P(f)  =  4p'  -  3p'  . 
b.  Five- jug  line  by  Equation  D5. 

i+1  /  5  -  i\  1 

,(f,-£<-,)  ( i  r  ■ 

recognize  that  |j?/2j  =  2,  i.e.,  retain  the  quotient. 


P(f)  =  (-D 


1  +1  /  5  -  1  \  1  2+1  (3  -  2  \  2 

(,)’■*  '-1’  (  JP'  ’ 


(DIO) 


(Dll) 


(D12) 


2 

P(f)  =  4p'  -  3p’  .  (DIO) 

Both  expansions  reduce  to  the  same  equation  for  the  probability  of  failure 
of  a  5- jug  line.  Evaluating  Equation  DIO  with  p1  =  (p/2)**2  =  (0.135/2)**2 
=  0.00456,  P(f)  =  0.01816  cs  0.0182  as  in  Table  1,  n  =  5. 


TABLE  D1.  JUGCOM  Program. 


PROGRAM  JUGCOM 

PROGRAM  TO  CALCULATE  the  PRO  wa?:  Il.it  OF  failure  ro 
D£7  UNATF  DUE  T 0  C it N E  E P  fi 0 F  <  m  I ...  ~  L  A 7  ED  JUG  S  .• 

IMPLICIT  DOUBLE  PRE  C I S  I C-N  •  A-h  0 -2  i 
WP I TE • 6,ii 

1  FORMAT  1  /  '  CALCULATION  Oc  PROBABILITY  OF  FAILURE  TO  DETONATE  •' 
i • DUE  TO  CRANE  ERROR  (MISPLACED  JUGS)'// 

* ' WHEN  PROMPTED,  ENTER  N.  THE  NUMBER  OF  JUGS  IN  A  ROW,/ 

+ 'AND  ® ,  THE  PROBABILITY  OF  AN  INDIVIDUAL  JUG  BEING  L A  CD'/ 

* ' OUT  SIDE  ITS  ERROR  CIRCLE.'// 

*  '  N  MUST  BE  A  POSITIVE  INTEGER  GREATER  THAN  OP  EQUAL  TO  £  .  ' / 

* ' AND  P  MUST  BE  GREATER  THAN  OR  EQUAL  TO  0  AND'/ 

* 'LESS  THAN  OR  EQUAL  TO  1  .  ' // 

*  ’  ENTER  ZEROES  FOR  N  AND  P  WHEN  NO  MOPE  DATA  ARE  TO  BE  ENTERED 
£  READ(5,*)N,P 

I F  (  N  EQ  .  0  )  GO  TO  5 
P£-  I  P  /£  .  0  2 

NTERMS=N/£ 

PTERM=~ 1 . 0 
PRFAIL-0  0 
DO  3  1=1,  NTEP.MS 

PTERM  =  -PTERM*P2*  (  (  NU  0 -?  .  0  *  I  )  /  (  N+ 1  0-1  >  )  =M  <  N  +  2 . 0-£  .  0*  I)  /  I  ) 

3  PPFAIL=PRFAIL+PTERM 
WR I TE ( 6 ■ 4 ) N , P , PRFA I L 

4  FORMAT (/'N  =' ,T16, 14/ 

*'P  *' ,T14,F6  4// 

* ' PROBABILITY ' / 

* 'OF  FAILURE  =  ' , F& . 4/) 

GO  TO  2 

5  STOP 
END 
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D.2  Exact  calculation,  JUGREC. 


a.  Union  of  Sets  A  and  B. 


b.  Intersections  of  Sets  A  and  B 


c.  A  and  B  as  Mutually  Exclusive  Sets. 


Figure  D-1.  Venn  Diagrams  for  Any  Two  Sets. 

Figure  D-1  shows  Venn  diagrams  of  any  two  sets  A  and  B.  In  Figure  D-la 
the  sets  overlap  because  they  have  some  elements  in  common  and  in  D-lc  they 
have  no  common  elements.  Every  element  of  the  black  bordered  "universe"  is 
either  in  "A  or  B",  written  A  U  B,  or  is  in  "not  (A  or  B)",  written  A  U  B. 
In  Figure  D-lb,  the  sub-regions  of  D-la  are_identif ied  as  elements  in  one 
set  and  not  in  the  other  set,  A  0  B  or  B  0  A;  the  overlap  region  of  A  and  B 
is  identified  as  A  O  B. 


By  writing  equivalent  expressions  for  the  region  of  A  U  B,  the  very  same 
elements  are  contained,  which  makes  the  sets  equal.  Inspection  of  Figure  D-lb 
and  D-lc  shows  that  the  region  identified  as  A  U  B  in  D-la  can  be  written  as 


A  U  B  =  A  +  B  -  (A  O  B), 


(D13) 


A  U  B  =  A  +  (B  H  A). 


(D14) 
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In  JUGFAE,  a  2-jug  line  involves  just  the  set  ,  so  Equation  D1  gives  failure 
for  n  =  2,  P(f)  =  P(E.,)  =  p*. 

For  n  =  3,  P(f)  =  P^ijEg)  =  P(E1  +  [E^E.,]),  (D15) 

=  pCe^  +  p(e2h¥1)> 

=  P(E.,)  +  p(e2)  =  2p '  , 

where  Equation  D14  and  the  additive  probability  of  mutually  exclusive  events 
and  the  fact  that  an  element  in  "not  E^"  has  to  be  in  E2  were  used. 

With  Equation  D15  as  a  model,  the  extension  to  a  4-jug  line  is  done  if 
in  place  of  E^  we  write  E^UE2>  and  in  place  of  E2  we  write  E^. 

For  n  =  4,  P(f)  =  P(E1UE2UE3), 

=  P(E1UE2)  +  P(E3n[E1uE2]).  (D16) 

All  manipulations  are  done  on  the  sets  in  the- second  term  so  we  work 
with  just  that  much  then  join  the  intermediate  result  to  Equation  D16. 

E  n[E  UE  ]  =  E  n[E  HE  ],  De  Morgan's  theorem, 

3  12  3  12 

=  E  n  E  flE  , 

3  1  2 

=  EHEnEflE  , 

3  13  2 

=  [E  HE  ]  n  [E  n  E  ], 

3  1  3  2 

=  [e^e.,]  n  e3, 

=  E3nE1HE3, 

E3n[EjuEj]  =  E3nfr 

So  Equation  D16  becomes 

p(f)  =  p(e1(je2)  +  p(e3Oe1),  (D17) 

E2)  +  P(E3)P(E1 ),  independent  events 

P(E^UE2)  +  P(E3)[l  -  P(E^)],  complementarity 

P(fail  w/4  jugs)  =  P ( f a i 1  w/3  jugs)  +  p'[l  -  P(fail  w/2  jugs)].  (D18) 

Equation  DIB  was  explicitly  derived  because  a  4- jug  line  is  the  shortest 
line  for  which  the  failure  equation  is  seen  as  a  specialization  of  the  general 
form.  We  go  on  to  derive  the  general  form,  but  the  symbology  is  tedious  to 
follow.  The  reader  may  prefer  to  note  that  the  result,  Equation  D25, 
subsumes  Equation  D18. 


A  D  A  =  A, 

association, 

E3  and  E2  do  not  intersect. 
e3  is  a  subset  of  "not  E./'. 

drop  braces, 

per  third  step  above. 
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Recall  from  Equation  D1  that  the  probability  of  failure  P(f)  for  n  jugs  is 

P(f )  =  P(E  UE  U  ...  U  E  UE  ), 

1  2  n-2  n-1 

=  P[(E  UE  U.--UE  UE  )UE  ],  association 
1  2  n-3  n-2  n-1 

=  P{ (E  U  E  U...UE  U  E  ) 

1  2  n-3  n-2 

+  [E  (E  U  E  U  •  •  •  U  E  UE  )]},  by  Equation  D14- 
n-1  1  2  n-3  n-2 

The  whole  expression  inside  braces  should  be  thought  of  as  two  sets, 

E  . ,  and  a  set  which_is  all  the  united  sets  inside  the  parentheses; 
its  form  is  {A  +  B  O  A}  which  is  the  sum  of  two  sets  that  are  always 
mutually  exclusive  (recall  Figure  D-1 ) .  So  the  probability  of  occurrence 
is  the  sum  of  the  separate  probabilities  of  occurrence,  as  was  already  used 
at  Equation  D16.  Then  we  can  write: 

P(f)  =  P(E  UE  U-..  UE  UE  ) 

1  2  n-3  n-2 

+  P[E  n  (E  UE  u...  U  E  u  e"  )].  (D19) 

n-1  1  2  n-3  n-2 

All  manipulations  are  performed  on  the  sets  in  the  second  term,  so  we  work 
with  it,  analogously  to  the  three-event  manipulations  after  Equation  D16, 
then  join  the  intermediate  result  to  Equation  D19: 

E  Pi  [  (E  UE  U...  U  E  UE  )]  =[E  H  (E  O  E  H  . . . 

n-1  1  2  n-3  n-2  n-1  12  (D20) 

Of  HE  )  ] ,  De  Morgan 
n-3  n-2 

=  E  OEHE  n...  nl  HE  ,  drop  parentheses 

n-1  1  2  n-3  n-2 

=  E  OE  OE  nEn...OE  HE  HE  HE  ,  AOB  HC  = 

n-1  1  n-1  2  n-1  n-3  n-1  n-2  A  H  B  H  A  H  C 
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since 


=  (E 

n 

e  )  n 

(E 

HE  )  n  ... 

n  (e  n 

E 

) 

n (E  HE  ) 

n 

-i 

i 

n- 

-1  2 

n-1 

n- 

-3 

n-1  n-2 

=  (E 

n 

e  )  n 

(E 

he  )  n 

n  (e  n 

E 

)n  E  , 

n- 

i 

n- 

-1  2 

n-1 

n- 

-3 

n-1 

the  last 

parentheses  are  dropped  because 

E 

n- 

-1 

^  ^n-2  =  ^n-1 

and  En 

_2  do  not 

intersect, 

=  (E 

n 

e  )  n 

(E 

he  )  n 

n(E  n 

E 

) 

9 

n^ 

-i 

1 

n- 

-1  2 

n-1 

n- 

-3 

where  the  last  Sn_^  is  dropped  because  it  is  already  intersected. 

For  the  same  reason  as  previously  stated,  all  but  one  En_.j  can  be  dropped. 

=  e  n e  nl  n...  ns  , 

n-1  1  2  n-3 

=  E  n  (E  HE  n  ...  HE  )  ,  association 
n-1  1  2  n-3 


=  E  n  (E  UE  U  ...  u  E  ),  De  Morgan 
n-1  1  2  n-3 

So  Equation  D20  becomes, 

(E  Pi[(E  UEU...  U  E  uF  )]  =  E  n  (E  UE  U-  UE  )•  (D21 ) 

.n-1  12  n-3  n-2  n-1  1  2  n-3 

Equation  D21  is  substituted  into  Equation  D19  and  it  becomes, 

P(f)  =  P(E  U  E  U...  U  E  UE  )  +  P(E  O  [  (E  u  E  U  •  •  •  U  E  )]•  (D22) 

1  2  n-3  n-2  n-1  1  2  n-3 

The  sets  in  the  second  term  are  independent  and  since  P(A  ft  B)  =  P(A)P(B) 
when  A  and  B  are  independent, 

P(f)  =  P(E  UE  U  •  ♦  •  U  E  UE  )  +  P(E  )[P(E  UE  U  ...UE  )],  (D23) 

1  2  n-3  n-2  n-1  1  2  n-3 

by  complementarity, 

P(f)  =  P(E  UE  U...  UE  UE  )  +  P(E  )[(l  -  P(E  UE  U...  UE  )],  (D24) 
1  2  n-3  n-2  n-1  1  2  n-3 

Equation  D24  is  analogous  to  the  equation  before  (D18)  and  with  that  one  and 
Equation  D1 ,  we  can  write, 

P(fail  w/n  jugs)  =  P(fail  w/n-1  jugs)  +  ( p ' ) [ 1  -  P(fail  w/n-2  jugs)],  (D25) 
2 

where  p'  =  (p/2)  ,  p  =  probability  of  one  jug  being  mislaid. 

Equation  D25  is  the  basis  of  the  program  JUGREC. 
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o  c_* 


TABLE  D2.  JUGREC  Program. 


PROGRAM  JUGREC 

PROGRAM  TO  CALCULATE  THE  PROBABILITY  OF  FAILURE  TO 

DETONATE  DUE  TO  CRANE  ERROR  (MISPLACED  JUGS) 

IMPLICIT  DOUBLE  PKEC I S I ON < A-H , O-Z » 

WRITE! 6, 1 ) 

1  FORMAT (/'CALCULATION  OF  PROBABILITY  OF  FAILURE  TO  DETONATE'/ 

* ' DUE  TO  CRANE  ERROR  (MISPLACED  JUGS)'// 

* ' WHEN  PROMPTED,  ENTER  N,  THE  NUMBER  OF  JUGS  IN  A  ROW,'/ 

♦•'AND  P,  THE  PROBABILITY  OF  AN  INDIVIDUAL  JUG  BEING  LAID'/ 
♦'OUTSIDE  ITS  ERROR  CIRCLE.'// 

*'N  MUST  BE  A  POSITIVE  INTEGER  GREATER  THAN  OR  EQUAL  TO  2,’/ 

♦'AND  P  MUST  BE  GREATER  THAN  OR  EQUAL  TO  0  AND'/ 

♦'LESS  THAN  OR  EQUAL  TO  1.'// 

♦'ENTER  ZEROES  FOR  N  AND  P  WHEN  NO  MORE  DATA  ARE  TO  BE  ENTERED.'/) 

2  R  E  A  D  (  5  ,  *  )  N ,  P 

I F ( N . EG . 0 ) GO  TO  5 
P£  = (P/2 ) ♦♦£ 

PRF  A I L£=  0 
PRFAILt -0 
PRFA I L=0 
DO  3  1=2, N 
PRFAIL2=PRFAIL1 
PRFAIL1 =PRFAIL 

3  PRFAIL^PRFaILH-PSM  1  -PRFAIL2  ) 

URITE(6,4JN,P,PRFAIL 

4  FORMAT ( / ' N  = ' , T1 6 , 14/ 

♦ ' P  = ' , T1 4 , F6 . 4// 

♦ 'PROBABILITY'/ 

♦ ' OF  FAILURE  =  ' , FA  4/ ) 

GO  TO  £ 

5  STOP 
END 
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^•3  Approximate  calculation.  A  quick  and  easy  approximation  can  be 
ased  if  Proceed  as  lr  all  of  the  events  E1 ,  E2 . En-1  were  independent. 

As  before  in  Equation  D1 ,  for  an  n-jug  line  the  probability  of  failure  is 


P(f)  =  P(E  UE  U...  U  E  ), 

1  2  n-1 

=  1  -  P(E  UE  U  . . .  u  E  ), 

1  2  n-1 

=  1  -  P(E  nl  fl...  HE  ), 

1  2  n-1 

-  1  -  P(E  )P(E  )  ...  P(E  ), 

1  2  n-1 

-  1  -  { [ 1  — P ( E  ) ] C 1  — P ( E  )]  ... 

1  2 

[l-P(E  )}, 
n-1 

n-1  2 

-  1  -  ( 1 -P ' )  ,  p'  =  (p/2) 

2  n-1 

P(f)  a  1  -  [1  -  (p/2)  ], 


complementarity 

(D26) 

De  Morgan's  theorems 

(D27) 

incorrectly  assumes  all 
events  are  independent 

(D28) 

complementarity 

(D29) 

(D30) 

(D31) 

where  p  is  the  probability  of  mislaying  one  jug. 


This  approximation  tends  to  underestimate  the  probability  of  failure. 
For  large  values  of  p  (close  to  one)  the  error  is  substantial  for  some 
values  of  n  but,  for  the  value  of  p  considered  in  this  work  (p  =  0.135), 
the  error  will  not  exceed  0.001  for  any  value  of  n.  High  values  of  p  are 
not  plausible  anyway,  as  JUGFAE  could  not  succeed.  The  approximation  in 
Equation  D31  would  be  suitable  if  some  quick  calculations  were  needed  in 
the  field,  where  sophisticated  computing  facilities  are  not  available. 
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